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Abstract 

In this paper, we present an algebraic formalism inspired by Butcher's B-series in numerical analysis and 
the Connes-Kreimer approach to perturbative renormalization. We first define power series of non linear 
operators and propose several applications, among which the perturbative solution of a fixed point equation 
using the non linear geometric series. Then, following Polchinski, we show how perturbative renormalization 
works for a non linear perturbation of a linear differential equation that governs the flow of effective actions. 
Furthermore, we define a general Hopf algebra of Feynman diagrams adapted to iterations of background 
field effective action computations. As a simple combinatorial illustration, we show how these techniques 
can be used to recover the universality of the Tutte polynomial and its relation to the g-state Potts model. 
As a more sophisticated example, we use ordered diagrams with decorations and external structures to solve 
the Polchinski's exact renormalization group equation. Finally, we work out an analogous construction for 
the Schwinger-Dyson equations, which yields a bijection between planar <fi 3 diagrams and a certain class of 
decorated rooted trees. 
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1 Introduction 

The last century has witnessed the development of Quantum Field Theory (QFT) as the adequate framework 
for the formulation of many physical theories, ranging from elementary particles to condensed matter physics. 
Historically, it was first devised to describe the interaction of charged particles with the electromagnetic 
field, leading to Quantum Electrodynamics (QED). However, QED was shown to be plagued by a surprising 
disease: When expanded in powers of the electric charge e, corrections to the leading terms turn out to be 
divergent. The occurrence of these infinities can be traced to virtual processes involving very high energies. 
Then, if we limit the energy these processes involve to a scale A and choose the electric charge e(A) as a 
suitable function of A, all the divergent quantities can be made finite, provided they are expressed in terms 
of a few low energy parameters to be determined by experiment, like the masses of the particles and the 
strength of their interactions. This is the basic idea of renormalization: Processes occurring at a given scale 
must be formulated using parameters adapted to these scale. 

Later on, QFT progressively invaded almost all fields of theoretical physics. Renormalizability was a 
guideline for the construction of the standard model of elementary particles while critical phenomena in 
condensed matter were also understood using renormalization. More in general, QFT describes fluctuating 
systems with continuous degrees of freedom and is best understood through paths integrals of the type 

f[Dcf>]e- s M, (1) 

where the integration is over a space of functions. When considered as a perturbation of a Gaufiian integral, 
this path integral is expanded over Feynman diagrams. Because of the continuous nature of the path integral, 
some of these diagrams yield divergent quantities, as dictated by dimensional analysis, and can only acquire 
a meaning through renormalization. In a certain sense, it can be said that the founding fathers of QFT 
(Dirac, Feynman, Schwinger, Tomonoga, Dyson, Weinberg, t'Hooft, Veltman, Wilson and many others) 
have achieved an infinite dimensional analogue of Newton's work: Driven by fundamental physics, they 
devised a new form of calculus based on path integrals, Feynman diagrams and renormalization. 

However, it is fair to say that this new calculus is not yet fully understood. In particular, it hides an 
algebraic structure analogous to diffcomorphisms, as uncovered by Connes and Kreimer a decade ago (see 
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P], [2] and [3] as well as the first chapter of the book [4] and [5] for recent reviews). To be specific, divergent 
Feynman diagrams generate a commutative Hopf algebra whose associated group is a refinement of the group 
of diffeomorphisms of the coupling constant. This formalism is versatile enough to encompass the infinite 
dimensional analogue of the change of parameters provided by the computation of the Wilsonian effective 
action, as we show in this paper. Our contribution is organized as follows. 

In the first section, we give a terse account of renormalization. We first recall the elementary particle 
physics point of view and introduce Feynman diagrams. Then, we introduce Wilson's viewpoint, first for spin 
systems and then for general QFT. This is textbook physics and can be found in modern texts on QFT (to 
quote only a few authors who strongly influenced our presentation, see the textbook [B] by Zee for an eclectic 
overview of QFT applied to particle and condensed matter physics, the treatise by Zinn- Justin [7] for a 
general treatment of QFT in the path integral formalism and the review by Rivasseau [8] for an introduction 
to renormalization, accessible to a large audience). We close this section by recalling a renormalization group 
proof of the central limit theorem of probability theory that may help mathematicians grasp Wilson's ideas. 

The next section is devoted to rooted trees and their Hopf algebra. These were introduced as funda- 
mental tools in the numerical analysis of a differential equation by Butcher |5J. We revisit this construction 
and interpret trees as indices for power series of non linear operators, with composition law given by the 
convolution product. We illustrate this construction with the non linear geometric series that provides a 
pertubative solution of a general fixed point equation. This can be used to resum tree like structures as 
shown on a elementary example based on Gaufiian random matrices. 

In the last section, we deal with renormalization. We first formulate perturbative renormalization for a 
general renormalization group differential equation. Then, we introduce a general Hopf algebra of Feynman 
diagrams that encodes successive effective action computations. We illustrate the use of this Hopf algebra 
on two examples. First, we give a new proof of the universality of Tutte polyomial and of its relation to 
the Potts model. Then, we use the Hopf algebra to solve Polchinski's equation that describes the flow of 
effective actions in QFT. Finally, in the last section we investigate an analogous structure based on the 
Schwinger-Dyson equations. As a byproduct, we obtain a new bijection between planar (j) 3 diagrams and a 
certain class of decorated rooted trees. 



2 Basics of wilsonian renormalization 

2.1 What is renormalization? 

In its most general acceptance, renormalization can be defined as the general methods that implement the 
change of parameters in the description of a physical system necessary in order to take into account the 
indirect effects of unobserved degrees of freedom, usually living on a much shorter length scale. The simplest 
example we can think of is the propagation of light in a medium, as described by geometrical optics. Whereas 
in the vacuum the speed of light is c, in some medium it is reduced to c/n, with n the refractive index, as a 
consequence of the interaction of light with the atoms of the medium. While the latter are not described by 
geometrical optics, their effect on the propagation of light is simply taken into account by the substitution 
c — > c/n. 

As a second example, let us consider the motion of a ping-pong ball of inertial mass mi nort , gravitational 
mass m grav and volume V in a perfect fluid of density p, as popularized by Connes jj. A straightforward 
application of Newton's law of dynamics yields its acceleration through 

m ine rta = m grav g - pVg, (2) 

where the first term is the ball's weight and the second one Archimedes' buoyant force. If the ball is ten 
time lighter than the equivalent volume of water, then it experiences an acceleration a = —9 g, which is 
obviously nonsensical. The hole in the previous reasoning is that we did not take into account the inertia 
of the surrounding fluid. Using a suitable ansatz for the fluid's velocity field, it can be shown jl] that the 
latter implies a shift of the inertial mass of the ball 

TOincrt — > minert + — , (3) 

which yields the more acceptable acceleration of a = — 3/2 g. Both these elementary examples illustrate 
the main feature of renormalization: In a simple theory (geometrical optics, Newton's law of dynamics), 
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the effect of unobserved degrees of freedom, which should be described by a more sophisticated theory 
(electrodynamics, fluid mechanics), are taken into account by a mere change of parameters in the simpler 
theory. 

The very same type of ideas show off in QFT, which grew up from the marriage between special relativity 
and quantum mechanics. First, recall that one of the main lessons from special relativity is the equivalence 
between matter and energy, as neatly summarized by Einstein's celebrated relation E — mc 2 . This is 
currently observed in high energy physics experiments: Incoming particles (say leptons, like a pair e + e~) 
with a total given energy collide and result in some outgoing particles whose nature may be totally different 
(for instance, a pair qq). Of course, total energy is always conserved so that heavy particles may appear 
only if the incoming energy is high enough. 

On the other hand, quantum mechanics provides us with a description of the microscopic world. Its 
predictions, carefully verified by experiments, yield the probability of finding the system in a final state |/) 
knowing that it was initially in the state \i). The transition amplitude from an initial state to the final one 
is often expressed as a sum over all (unobserved) intermediate states \ri), 

(f\i) = J2(f\n){n\i). (4) 

n 

In QFT, these initial and final states represent a fixed number of elementary particles of given energy but the 
summation over all the intermediate states involves an arbitrary number of virtual particles of arbitrarily high 
energies. It is fundamental to note that none of the intermediate states are observed since the experiment 
only answers the question "What is the probability of finding the system in the state \f) at the end of the 
experiment if we know that it was initially in the state \i) ?". 

In most cases of interest, the sum over intermediate states diverges as the upper limit of the energy A of 
the virtual particles go to infinity. However, for the so-called renormalizable theories, it turns out that these 
infinities cancel when the parameters of the theory (masses m(A) and coupling constants 5(A)) are suitably 
chosen as functions of A, which may then be taken to infinity leaving physically relevant quantities finite. 
This is the process of renormalization, as it occurs in high energy physics: The effect of unobserved virtual 
states of high energy is taken into account in some of the parameters of the theory. In a certain sense, the 
cut-off A sets a limit for the validity of our field theory and the sole effect of the very high energy modes 
beyond A is to adapt the parameters of the theory to energy of the process we consider. Alternatively, this 
can also be considered as an ignorance of the very small scale structure of space-time since high energy 
correspond to small distances by Fourier transform. 



2.2 Feynman diagrams and perturbative renomalization 

As far as renormalization is concerned, QFT is best studied in the framework of the Euclidian path integral. 
In the simplest case, the fields are defined as functions <fi from space-time M. D to K and the dynamics is 
governed by the action functional, 

S[<j>] = J d D x (^-d^d^+\m 2 4 ) 2 + ^ N ^j , (5) 

where to is a mass and g a coupling constant for a A^-particlc interaction. At the quantum level, all the 
information is encoded in the Green's functions, 

/_ gM 
[D(j))e *■ (j){xi)---4>{x n ), (6) 

where the integration is over the space of all fields and Af is a normalization constant to be defined later. 
h — 6.02 1 -34 J-s is Planck's constant and measures the deviation of the quantum theory from the classical 
one. In most of the paper, we choose unit such that h = 1. Note that is is only the simplest model of a 
scalar field theory. More complicated action functionals are required to account for the real world physics: 
Spinors ip for fermionic particles and gauge connections A for their interactions, as is the case for QED and 
the Standard Model of elementary particles. Nevertheless, we restrict our attention in the sequel to the 
simplest example of a scalar field theory. 

To give a precise meaning to |6]), it is convenient to expand G(x±, . . . , x n , m, g) as a power series in g using 
Feynman diagrams. The latter are best introduced on simpler analogue, with the space of fields replaced by 
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a finite dimensional vector space. Thus, the equivalent of ^ is 

G il ,..., in (A,V)=Jf I tye-SMfa'-'fa, (7) 



with the action 

S(4>) = ^-A- 1 - ( t> + V( ( p). (8) 
The quadratic term is defined by a positive definite symmetric matrix A, 

t-A- 1 -cf> = Y f (A- 1 )i j cl> j( t>i, (9) 
and the interaction potential V is a polynomial in all the fields 

^) = E E 9jl ff L ^---4>i N - (io) 

N ii,...,ijv 

In this case we choose N^ 1 = det A/2tt and expand e~ v ^ as a power series in the couplings gi lt ...,i N . Thus, 
([7]) amounts to the computation of the average of a monomial using a GauBian weight and is given by Wick's 
theorem: It is a sum over all possible pairings of the variables in the monomial, each pairing of fa with <j)j 
being weighted by Ay. For example, 

fa<t>j4>kfa — * AijA k i + A ik Aji + AuAjk. (11) 

Then, each term of the expansion of ^ is associated with a diagram with n external legs and vertices of 
valence N, 

_ G, 7 . (A, V) 
G iu ... tin (A,V) = ]T U -'o ■ (12) 

7 ' 

The contribution of each diagram is computed using the Feynman rules: 

• associate the indices i\, . . . ,i n to the external legs and indices jk to the internal half edges; 

• associate a matrix element Aj k j l to any edge connecting the indices jk and ji\ 

• associate a coupling —gj 1} ...,j N to a iV-valent vertex whose half edges have indices j\, . . . , j^\ 

• sum over all the indices jk- 

Besides, on has to divide by the symmetry factor S 7 which is the cardinal of the automorphism group of the 
diagram, leaving the external legs fixed. For example, 



r, Ai 1 j 1 gj 1 j 2 .j 3 Aj 2 j i Aj 3 j 5 gj 4 j 5 j 6 Aj 6t i 2 . (13) 
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31,32-33 
34*35,36 



This simple finite dimensional model already captures some important algebraic aspects of perturbation 
theory as will be discussed in the last part devoted to the Hopf algebras based on Feynman diagrams. 

At a formal level, the Green's functions ^ can be computed as a power series in g by replacing fa by a 
function x i— > fax), the matrix element Ay by the propagator 



d D p e ip-(x-y) 

(2ir) D p 2 + m 2 



K(x,y)= I t^-^t^. (14) 



and V(fa) by the interaction term 



-jL f d D x^{ X ) = j-( ^...fP^(2n) D S(p 1 + ---+p N )fap 1 )...fap N ), (15) 

with 

fap) = / d D X e~ ip x fax) (16) 
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the Fourier transform of 4>{x). Heuristically, the Fcynman diagrams can be thought of as quantum mechanical 
processes with particles on their external legs and virtual particles of momenta p propagating on the internal 
lines. It is important to notice that although momentum is conserved at each vertex and along each line, 
the particles that propagate along the loop^Jmay have arbitrary momenta. 

The integral along loop momenta extend on the whole of R"° and may led to divergencies since the Fourier 
transform of propagator does not decrease fast enough at large p. A connected diagram with / internal lines, 
V vertices and L = I — V + 1 loops behaves in dimension D like p u with u> = LD — 21 the superficial degree 
of divergence. Using NV = 21 + E, the latter can be written as 

V (ND „ .A „A D 



V\—--D-N\+E\l--\+D (17) 

When the coefficient of V is strictly positive, divergencies occur for any Green's function and such a theory 
cannot be renormalized. If this coeffcient is strictly negative, then there are only a finite number of divergent 
diagrams. These are the superrenormalizable theories which are nevertheless of limited interest in physics. 
Finally, the critical situation defines the renormalizable theories, like </> 4 in D = 4 or </> 3 in D = 6. In this 
case, the divergencies occur only for the diagrams with E < N and may absorbed into the coefficients of a 
polynomial interaction of degree N. At the notable exception of gravity, all the interactions of elementary 
particles can be formulated using renormalizable interactions. 

To prevent the propagation of Fourier modes of momenta > A, let us alter the propagator by introducing 
a cut-off A, 

r°o t d D p 



K(x,y) - K A (x,y)= \ da \ ^<-v) e -(i +»» ). (18) 



This procedure is known as the regularization and can be performed in various ways: Besides the method 
used here, one could also discretize the theory on a lattice or evaluate the diagrams in complex dimension 
z and recover the divergences as poles when z — > D. In principle all these methods are equivalent but we 
restrict ourselves to the momentum space cut-off presented here since it is suited to the Wilsonian point of 
view we adopt in this paper. 

For a renormalizable theory like the 4 theory, one can trade the parameters g and m for some cut-offs 
dependent ones go (A) and mo (A) and further introduce an additional wave function renormalization Z(A) in 
such a way that Z? (A)G\(xi, . . . ,x n ,mo(A),go(A)) admits a finite limit when A — > oo. To obtain definite 
physical predictions, the bare parameters go(A) and mo (A) and the wave function renormalization Z(A) must 
be determined in terms of normalization conditions involving renormalized parameters m r and g T measured 
at a low energy scale /i. Thus, we define the renormalized Green's functions as 

G T (xi,...,x n ,m r ,n,g T ) = lim Z% (A, m r , g T , fi)G A (xi, . . . , x n , m (A, m T ,g T , fi),g {A, m T ,g T , n)) (19) 

A— >oo 

Note that we are dealing here with perturbative renormalization only, so that the previous equality must be 
understood as an equality between formal power series in g r . In fact, <?o(A), mo (A) and Z(A) are themselves 
formal power series in g r that can be computed in terms of Feynman diagrams using the Bogoliubov-Parasiuk- 
Hepp-Zimermann (BPHZ) formula. Roughly speaking, the contribution of a divergent diagram with 2 or 4 
external legs to the renormalisation of the parameters is encoded in its counterterm C{^) which is determined 
recursively by the relation 

c(7) = t( E II^ttV)' (20) 

{t»,...} i iU h 

with T taking the divergent part of the diagram. This sum runs over all sets , including the empty one, 
of disjoint, divergent, one particle irreducible subdiagrams of 7 (i.e. diagrams that cannot be disconnected 
by cutting an arbitrary internal line). The reduced diagram on the RHS is obtained by shrinking each 7$ 
to a single vertex and finally take the divergent part of the whole sum, with a finite part determined by 
the normalization conditions at the low energy scale /i. In the framework of dimensional regularization, this 
operation is elegantly written as a Birkhoff decomposition for a loop in the space complex dimension with 
values in a group associated to a commutative Hopf algebra (see the work of Connes and Krcimer [5] and 
0)- 

4 What is called loop in the physics terminology adopted here is called a cycle by graph theorists. 
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Figure 1: One loop diagram for the four point function 




Figure 2: Chain of one loop diagrams for the four point function 



In the construction of the renormalized theory, the renormalizcd coupling constant g r and mass m r are 
the parameters of the theory which are directly accessible. They are determined by experiments performed at 
a low energy scale fi. However, if we consider the same theory but with renormalized parameters determined 
at a slightly different scale // , then we are performing a mere change in the parametrization of the same 
theory. Consequently, the renormalized parameters measured at different energy scale are not independent. 
For instance, the renormalized coupling constant depends on the scale through a differential equation of the 
type 

4,= m (21) 

The (3 function can be determined by requiring that the bare paramaters, which are the "true" parameters 
of the theory entering in the path integral, are independent of /i. Accordingly, the low energy renormalized 
theory has a subtle form of scale invariancc inherited from the behavior of the high energy virtual processes. 
The zeroes of (3 define the fixed points of the renormalization group and play a fundamental role. In 
particular, the origin g = is a fixed point and it is (UV) attractive if /3(g) < 0. This last case corresponds 
to the asymptotically free theories for which the coupling constant goes to at high energies. These theories 
are believed to be consistent non pertubatively, that is, the expansion of the Green's functions can be given 
a meaning beyond that of an asymptotic series at g = 0. In the other case, there are strong arguments 
against a non perturbative definition. 

For instance, at one loop level for the <p A model in dimension 4, the relation between the bare and the 
renormalized theory follows from the computation of the diagram given in figure [T] 

- 9r = -.90 + I3g 2 log A + O( 5o 2 ), (22) 
A* 

with f3 a positive constant. Here, g r is identified with the four point function evaluated on a specific 
configuration of momenta of order fj,, using the bare theory with a cut-off A. 

From this relation, one readily extracts the (3 function as (3(g) = fig 2 , so that the theory is not asymp- 



totically free. If we now attempt to get a non perturbative information from (22 1 by summing the geometric 
series associated to the chain of one loop diagrams depicted on figure [2] we get, as go > for the consistency 
of the path integral, 

9r = ;r < — ^r, (23) 

yr l + /3 5o log^ " /31ogf 

so that the g r — ► as A — > oo. This is the triviality problem of tfi 4 theory. Note that the argument would 
not apply if (3 were negative. 

To close this preliminary section on renormalization, let us give a brief historical sketch to show how 
these ideas gradually appeared in physics. 

• Early computations in QED 

It was already noticed by the founding fathers of QED in the thirties that the first non trivial terms 
of the perturbative expansion are plagued by infinities. However, certain combination of divergent 
quantities were shown to be finite, thus leading to successfully tested predictions. For instance, the 
anomalous magnetic moment of the electron, which is nowadays one of the most precisely tested 
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prediction of physics. This set the stage for the general proof of the renormalizabilty of QED in the 
fifties by Dyson, Feynman, Schwinger and Tomonoga. 

• Perturbative QED and its /3 function 

For QED with massless fermions an altered form of scale invariance was discovered: If we measure 
the low energy electric charge e at various energies /i, we describe the same physics provided we allow 
for a /i dependent charge e M . QED has a positive (3 function so that the electric charge increases at 
high energy. This also indicates a breakdown of (perturbative) QED at very high energy, which is no 
surprise since QED is not a complete theory as it does not take into account nuclear forces. 

• BPHZ (Bogoliubov-Parasiuk-Hepp-Zimmermann) formula 

Systematic combinatorial rules for renormalization were found in the sixities allowing to express the 
counterterm of a diagrams in terms of counterterms of its subdiagrams and divergent integrals. This 
is the content of the BPHZ formula and solved the problem of overlapping divergences. 

• Electroweak interactions 

In the late sixties, Glashow, Salam and Weinberg proposed a model based on spontaneously broken 
non-abelian gauge theories combining electromagnetic and weak interactions. Note that the latter was 
previously described by the non renormalizable Fermi Lagrangian. This electroweak model was later 
shown to be renormalizable by 't Hooft and Veltman and the predicted vector bosons were observed at 
CERN in 1983. Spontaneous symmetry breaking yielding massive vector bosons is achieved through 
the introduction of the Higgs field, still to be experimentally discovered at the time of writing (mid 
2008). 

• Strong interactions 

In the seventies, non-abelian gauge theories with a limited number of fermions were shown to be 
asymptotically free, meaning that the scale dependent coupling constant is driven to zero at very high 
energy. This is in agreement with experiments performed at SLAC where quarks inside the hadrons 
appear to behave as free particles at high energies. Accordingly, the strong force is also described by 
a renormalizable field theory, Quantum ChromoDynamics (QCD). However, at low energy the QCD 
coupling constant increases thus impeding a perturbative treatment. 

• Critical phenomena 

Critical phenomena occur in second order phase transitions where the correlation length tends to 
diverge, so that many degrees of freedom interact. In the seventies, the seminal work of Wilson led to 
an understanding of critical phenomena, making use of techniques of renormalizable QFT. 

• Effective field theories 

Later on, the Wilsonian viewpoint on renormalization shed a new light on high energy phenomenology. 
Indeed, most of the models currently used are thought of as effective field theories, yielding accurate 
physical predictions within a certain range of energy. This is also the point of view adopted in string 
theory: Whereas the latter is valid at an energy of the order 10 19 GeV (Planck's energy) or beyond, 
the standard model coupled to gravity is understood as a low energy effective action valid at a few 
hundred GeV. This is just as Fermi's theory that can be understood as a low energy effective theory 
derived from the electroweak model. 

• Constructive QFT 

In fact, the renormalization program in QFT we previously outlined only makes sense at the pertubative 
level. This means that physical quantities are computed as asymptotic scries in a small coupling 
constant but nothing is said about the convergence of this expansion. Constructive QFT is a branch 
of mathematical physics aiming at defining the theory using general principles and the information 
encoded in its perturbative expansion (see [2]). Here too, the Wilsonian point of view led to substantial 
progress. 

• Stochastic dynamics 



8 



QFT appear to be a powerful tool each time a system is modeled by a partial differential equation 
involving a random driving force. For example, the random growth of a surface is described by 
the Kardar-Parisi-Zhang equation whose solution is expressed as a path integral (see [5] for a nice 
pedagogical introduction). Then, renormalization theory provides some efficient tools to investigate 
the universal behavior of such a system. 

The ubiquity of renormalization theory in the previous situations can be understood as follows. In each 
case, the physical description is supposed to be valid on a wide range of length scales, including extremely 
small ones that are not directly probed by experiments. Then, renormalization is the set of methods that 
prescribes how the parameters of the theory have to be adapted to the degrees of freedom that live on a 
given scale to take into account the effect of the unobserved degrees of freedom that live on smaller scale (or 
higher energy). 



2.3 Coarse graining for spin systems 

The idea that renormalization is an adaptation of the theory to its natural scale mostly follows from the 
work of Wilson on critical phemomena. To illustrate how the Wilsonian point of view appeared in the study 
of critical phenomena in statistical mechanics, let us consider a system of spins cjj located on the sites of 
a lattice of spacing a. To simplify the discussion much as possible, let us assume that the spins can only 
take q different values on each site with a Hamiltonian H involving short range interactions (say a nearest 
neighbor one) that favors spin alignment. For instance, one can take 

ff(a) = -J^<5 CTi)CT ., (24) 

where the sum runs over all pairs of nearest neighbors and J is a positive constant. This is the celebrated 
Ising model for q — 2 and the g-state Potts model in the general case. 
The basic object of interest is the partition function 



e-^\ (25) 



where the sum runs over all spin configurations, (3 — ^ , with k — 1.381 1CP 23 J-K -1 Boltzmann's constant 
and T the temperature. For these models, there is a competition between the tendency of the spin to align 
in a way to minimize their energy and thermal fluctuations which tend to favor random orientations. Some 
specific systems like the two dimensional Ising model exhibit a phase transition with a high temperature 
T > T c disordered phase and a low temperature T < T c ordered one, with T c the critical temperature. As the 
temperature approaches the critical one, the correlation length tend to diverge, thus leading to long range 
collective phenomena. Moreover, its has been observed that critical systems fall into universality classes 
made of systems that have a common asymptotic behavior. For instance, in a given universality class, the 
correlation length may diverge as T — > T c with the same exponent, 

f ~iriW- (26) 



Wilson's renormalization group idea amounts to perform the summation in (251 step by step in order to 
end up with a much simpler description. Specifically, one successively averages over spin configurations on 
lattices of increasing spacing (see figure [3]) 



g o 1 (27) 



together with a renormalization group transformation of the Hamiltonian 



The latter is defined by 



H ^ H'. (28) 



e -m'(<r') = J2 e~^\ (29) 

T(a)=a> 
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Figure 3: Coarse graining 



in such a way that the partition function remains invariant 



-PH'(a') 



(30) 



As we perform the coarse graining transformation from a to a', the lattice size is doubled and it is 
convenient to also resize the unit of length by a factor of s = 2, so as to stay on the same lattice. Accordingly, 
the correlation length is divided by s and after sufficiently many iterations of the renormalization group 
transformation, one can hope to end up with a tractable problem with a correlation length of the order of 
the lattice spacing. Moreover, the sequence of Hamiltonians defined by 



undergoes some drastic simplifications if it converges towards a fixed point H* of 1Z. Let us further assume 
that the linearization of 1Z in the vicinity of H* can be diagonalized with eigenvalues that are written in the 
form Xi = s Vi and eigenvectors Oi. Starting with a Hamiltonian close to the fixed point, 



the iterations of 1Z lead us to distinguish to three cases. 

• Relevant terms associated with > tend to grow and one has to fine tune the initial Hamiltonian 
very close to the critical surface (i.e. the basin of attraction of the fixed point) in order to avoid 
divergences. Relevant terms govern the deviation from criticality and there is usually one of them 
proportional to T — T c . 

• Irrelevant terms corresponding to yi < die off very quickly and do not play any significant role. 
Information pertaining to them is erased so that universality is recovered. 

• Marginal terms with zero eigenvalue evolve very slowly but after a deeper study they finally end up in 
one of the two previous classes, apart from some exceptional cases. 

As an elementary illustration, let us determine the critical behavior of the correlation length as an inverse 
power of \T — T c \. We start with a system close to criticality, so that its relevant parameter t = T — T c with 
eigenvalue y is very small and its correlation length £ large. We now iterate n times the renormalization 
group transformation in such a way that the relevant parameter t' becomes of order unity and the correlation 
length £' close to the lattice spacing, 



H n+ i — TZ(H n ) 



(31) 




(32) 



a „ g „ A_ and i ^ t > „ g ny t 



(33) 



Eliminating s n we obtain the critical behavior of the correlation, 



1 



(34) 



\T-T c \v 



where we have discarded a slowly varying multiplicative constant . 
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To end this very brief introduction to Wilson's ideas, let us sketch a renormalization group proof of the 
central limit theorem in probability theory, which has the merit of explaining the universal nature of the 
GauBian probability law (see and [H]). Consider N — 2™ independent identically distributed (IID) 
random variables £j with probability density po, assuming that its variance is 1 and a vanishing expectation 
value. Recall that the central limit theorem states that 



c 



N 



N 



(35) 



converges in law towards the normal law with density 2 as N — ► 00. 

In the spirit of coarse graining, let us compute (jv by iterating n times the transformation that associates 
to two IID random variables £' and £" the new one defined by 



V2 ' 



The probability density IZ(p) of £ is expressed in terms of that of £' and £" as 

[ft(p)] (x) = V2 [ dyp(y)p(V2x - y), 



so that the probability density of £jv can be written as 

p„ = n n ( Po ). 

The renormalization map TZ admits the GauBian 



p*(x) 



-x 2 /2 



/2lT 



(36) 



(37) 



(38) 



(39) 



as a fixed point. To investigate the stability of this fixed point, let us linearize TZ in the neighborhood of the 
GauBian. Writing the perturbation as 



we have, up to order e 2 terms, 
K[p](x) = 



J_ -x 2 /2 



1+^ / dye-«\(y + x/V2) 



+ 0(e 2 ). 



(40) 



(41) 



The eigenvectors of the linearization of 1Z are obtained by taking e proportional to a Hermite polyomial H n 
This is easy to check using their generating function 



(42) 



The associated eigenvalues are A„ = 2 1_ 5, so that all the terms involving Hermite polynomials of order 
n > 2 are irrelevant. For n — 0,1,2, the coefficients of the expansion of e have to vanish in order to make 
sure that the iteration of 1Z takes the probability density towards the fixed point. They simply state that 
the probability distribution has to be normalized (n — 0), with vanishing expectation value (n = 1) and 
with variance 1 (n = 2). In fact, only n = and n = 1 is are relevant, n = 2 is marginal. Imposing 
these conditions on the probability density we start with is analogous to the fixing the otherwise divergent 
quantities in QFT by low energy measurements. 

2.4 Wilsonian renormalization in QFT 

The very same ideas can be applied to QFT, with the partition function replaced by the path integral 



(43) 
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The integration measure involves an UV cut-off A, restricting the integration over fields whose Fourier modes 
vanish for momenta above A. It is analogous to a system defined on a lattice and reflects the fact that rapid 
variations of the fields on spatial scales below a ~ A do not make sense in this context. 

In analogy with spin blocking, one integrates over fast modes in order to obtain a equivalent field theory 
with a lower cut-off A' = — . To proceed, let us separate the field <j> into its fast component <p" with momenta 
between A' and A and its slow one cj>' with momenta below A'. By performing the integration over </>", we 
obtain the low energy effective action 

e" s ' M = J [^'V,Ae- s ^"l, (44) 

The partition function may be computed either with S or with S', 

Z = J [D<f>] A e~ s ^ = J [ D( j>>] A , e- s '^. (45) 

The equality between the two terms follows form a substitution of S' in terms of S followed by the change 
of variables <f> = <f>' + <p" . The remaining integration is cancelled by a normalization factor. These are self- 
understood constants in all formulae that do not depend on the actions and are defined by imposing that 
effective actions and partition functions are trivial for trivial actions. 

Besides, this remains true if we insert in the path integral any functional of <f>' , so that one can conclude 
that low energy physics remains unchanged if we simultaneously lower the cut-off and replace S by S' . 
All the low energy effects of the physics at the intermediate scale A' have been encoded in S' , which may 
therefore have a rather complicated form even if S is simple. 

As for spin systems, one is led to iterate the renormalization group transform to obtain a sequence of 
effective actions with decreasing cut-offs. In the context of QFT, it more transparent to adopt a continuous 
formulation, leading to a family of cut-off dependent effective actions S Al governed by a differential equation 

A^ = 0(A,5 A ). (46) 

The action with started with is encoded in the initial condition, S Aa = Sq, valid at a very high energy scale 
Aq. To derive the renormalization group equation (46), let us note that S A must be obtained from S Ao by 



integrating over all the modes between A and Ao. This is achieved in the path integral approach by 

e -S A [<t>] _ J ^£) x ^ e -\x-Al\ Q x-SA [x+4>] ^ ( 47 ) 

where Aa,a i s a smooth cut-off dependent propagator that implements a smooth analogue of the integration 
over the fields with Fourier modes between A and Ao. For example, we shall often use a cut-off inspired by 
Schwinger's representation 

Aa. Ao (p, q) = S(p,q)[™ da e" Q f 2 , (48) 

(A ) 2 

with the mass term relegated in the effective action. Note that A plays the role of an IR cut-off whereas Ao 
is an UV one. By differentiating with respect to A, we get Polchinski's equation (see 1 OJ ) , 

t dS A If,, K dA AAor , / S 2 S SS 8S \ 

A— -p- = - / dpdqk Jr: (p,q) — , (49) 

OA 2 J OA ^ ,HJ \S4>(p)84>(q) 8cb(p) 8<b(q) J V ' 

which is conveniently captured by Feynman diagrams as illustrated on figure [4] Although A\ Ao depends on 
both the UV and IR cut-offs, its derivative only depends on A. 

Here we have use the Fourier transform <p{p) of the field 4>{x) but an analogous equation can be written 
in position space, 

K dS A If,,, dA A Ao/ . ( S 2 S 8S 8S \ 

~dA =- 2 J^dyA — ( X ,y) [j^^ - W)JW) ) ■ (50) 



with 



A A ,Ao(:r,y)= / -^^ ^ A AM {p-p). (51) 
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Figure 4: Diagrammatic form of Polchinski's equation 



To close this terse section on renormalization, let us draw a few high energy physics consequences from 
the Wilsonian approach. In this context, Ao is the UV cut-off introduced in order to render the Feynman 
diagrams convergent and Sq identifies with the bare action. S\ is a low energy effective action, which 
has to be further used in Feynman diagrams having all their internal momenta below a fixed low energy 
scale A. As we try to get rid of Ao by taking it to infinity, some of the parameters entering in S\ tend 
to diverge. The latter have to be fixed by low energy measurement and renormalization simply amounts 
to choosing a bare action Sq in such a way that the flow will meet these low energy conditions. Since the 
limit ^ — > oo corresponds to an iteration of a large number of renormalization group transforms, this is 
very close to what has been done for critical phenomena. From this point of view, one can understand 
the ubiquity of renormalizable theories depending on a few paramaters as an analogue of universality in 
statistical mechanics. Whatever Sq we start with, provided it is natural (i.e. its parameters are numbers ~ 1 
times powers of Ao dictated by dimensional analysis) , we always end up with up with renormalizable theories 
at energies A< Ao. By the same token, universality implies that the renormalized theory is independent of 
way we implement the cut-off. 

Besides, this analogy can be pushed further: All the renomalizable interactions that appear in high 
energy physics can be thought of as low energy effective theories derived using a renormalization group 
flow from a yet unknown theory, valid at a very high energy scale Ao. The latter theory therefore plays in 
particle physics a role analogue to lattice models in condensed matter physics. As such, it may not even be 
a field theory and rely on some utterly new ideas, like string theory, spinfoam models or noncommutative 
geometry. The important point is that it comes equipped with a high energy cutoff Ao, of the order of 
Planck's length, which is physical and should not be taken to infinity. At energies slightly below Ao, the 
physics is well captured by a QFT whose field content, action and symmetries are postulated, but should 
determined by the unknown theory. Starting with this action, we move towards lower energies A using the 
renormalization group equation and end up with an effective action determined by a few relevant parameters 
to be experimentally determined. Irrelevant terms corresponding to non renormalizable Lagrangians are 
also present, though very small, and are innocuous because the cut-off remains finite. Therefore, there is no 
special problem with gravity, the latter being simply the leading irrelevant correction to the Standard Model, 
and the Einstein- Hilbcrt Lagrangian may be itself corrected by higher order terms. Of course this does by no 
means solve the problem of quantum gravity, the later presumably amounts to finding the unknown theory 
we should start with. 

One of the main lessons from the Wilsonian approach to renormalization in high energy physics is that 
the cut-off acquires a physical meaning: It sets the upper limit of validity of the theory. This is in sharp 
contrast with the point of view adopted a few decades ago where the cut-off was a mere scaffolding needed 
to build the theory and that has to be removed as soon as possible. 

This approach also raises another question: How do we make sure that the relevant paramaters, like 
the masses, are small and not of the order of Ao? For lattice systems, an operator can tune the bare 
Hamiltonian sufficiently close to the critical surface so that we end up with finite results even after a large 
number of iterations. In particle physics, finding such a process remains an open question and is particularly 
important for a scalar boson like the Higgs field. It seems that such a problem can be solved by introducing 
supersymmetric particles but at the tome of writing (mid 2008) no such particles have been found. 

3 Rooted trees and power series of non linear operators 
3.1 Rooted trees in numerical analysis 

Soon after Connes and Kreimer presented an algebraic framework base on rooted trees for renormalization 
PQ, it was realized by Brouder |14j that perturbative renormalization has some intriguing connection with 
the algebraic techniques used in the numerical analysis of differential equations. 
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Figure 5: A rooted tree with one branching 



as 



To begin with, let us expand in powers of s — so the solution of the differential equation in 

X(x), x(s ) = x , 
(s - s ) n d n x 



dx 
ds 



ds n 



s=s 



(52) 
(53) 



To proceed, we compute the successive derivatives of x with respect to s, 

X i 

(J.ri 



dx l 

ds. 
d 2 x % 

~df? 

d 3 x i 
Is* 

d A x i 
"ds^ 



dX i dX j 



^ dxi dx k 

j-k 

v dX l dXi dX k , 

l^^-^^x 1 



d 2 X l 



dxidx k 



X ] X k 



j,k,l 



dxi dx k dx 1 



d 3 X l 



dxi dx k dx l 



X j x k x i 



3 d 2 X* dX k xJxl 
dxidx k dx 1 

dx* x , xl 



dxi dx k dx 



(54) 



All the terms in this expansion turn out to be in one-to-one correspondence with rooted trees with at most 

.2 

four vertices. For example, the second term in the derivative %r§ corresponds to the tree depicted in figure 

El 

This is familiar from Runge-Kutta methods in numerical analysis. Indeed, a Runge-Kutta method is 
defined by a square matrix (o»j - )i<*j<n an d a vector (£»i) i<i<n of real numbers, arranged into a table 



(55) 



Cl 


an 


a\2 ■ 


ain 


C2 


0,21 


a-22 ■ 


■ a 2n 




a nl 


a n 2 ■ 


("Tin 




bi 


b 2 ■ 





where c t = J2i 



Then, the solution of the differential equation 

X(x), with x(sq) = xq, 



dx 
ds 



(56) 



is given at si = So + h by xi = xq + ^2i=i biX(yi), where y^ is determined by yi 



x 



When dij — for j > i, the computation of yi in terms of trees is straightfoward. Then, trees turn out to 



be useful to compare the result of Runge-Kutta method with the expansion (54 1. 

It is an amazing fact discovered by Butcher that Runge-Kutta methods can be composed One can 
perform a first approximate computation using a first method to obtain a: as a function of the initial condition 
xq. Then, we can consider i as a new initial condition and proceed to a different approximate computation 
resulting in x' . This amounts to a single computation of x' with initial condition xq, using a third method 
which is the product of the previous ones. 
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Figure 6: An admissible cut on the rooted tree of figure [5] 



3.2 Hopf algebra structure 

The occurrence of rooted trees in numerical analysis through Runge-Kutta methods and their composition 
law is conveniently described using the Hopf algebra of rooted trees. To begin with, recall that a tree is a 
connected diagram without loopj^jand a rooted tree is a tree with a distinguished vertex called the root. A 
rooted tree is oriented from the root to its terminal vertices, the leaves. Except in the very last paragraph 
devoted to a bijection between planar diagrams and trees, we always consider isomorphism classes of trees. 

Let TLt be the commutative algebra generated by all (isomorphism classes of) rooted trees. It admits a 
Hopf algebra structure defined on the generators as follows. 

• The coproduct A : TLt ~ > <8 TLt is 

A(i) =i(g>l + l<g)t+ P c (t) <E> R c (t) (57) 

c admissible cut 

A cut of an edge is admissible cut if any path from any leaf to the root is cut at most once. R c (t) is 
the connected component that contain the root after the cut, it is made of vertices located above the 
cut. P c (t) is the product of remaining trees corresponding to vertices below the cut with as new root 
those vertices immediately below the cut. For example, for the simplest tree with one branching, we 
have 



A (d^) =1 ®c/\ + c/\® 1 + 2#0 I + # * 0# - (58) 

• The counit e : TLt — > C is trivial, e(t) = unless t = 1. 

• The antipode S : TLt — * is defined as a sum over all cuts, 

S(t) = -J2(- 1 ) ncit)n o(t), (59) 

c cut 

with n c (t) the number of edges of the cut. 
Since the Hopf algebra TLt is commutative, its characters form a group Gt for the convolution product 

a */3=(a®/3)oA (60) 

with unit e and inverse a -1 = a o S. Recall that a character is an algebra morphism from TLt into a fixed 
commutative ring, 

J a(Xa + fib) = X a(a) + /j a(b) 

{ a(ab) = a(a)a(b), 1 > 

for any a,b € TLt and A,/ieC. Let us emphasize that the characters must take their values in a commutative 
ring, in order to check that that the product of characters is still a character. Most of the time we work 
with complex numbers but it may be useful to consider rings of functions. 

A character of TLt is entirely specified by its values on the trees, which may be arbitrary, so that it 
encodes an infinite sequence of complex numbers. In fact, Gt is a Lie group whose Lie algebra Qt consists 
in all infinitesimal characters. The latter are linear maps 8 : TLt —> C, satisfying the infinitesimal form of 
multiplicativity, 

5(ab) = e(a)6(b) + 6(a)e(b) (62) 



3 We use the QFT terminology: A loop is what graph theorists call a cycle. 
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for any a, b £ Ht- Alternatively, infinitesimal characters can be defined as derivations with values in the 
trivial bimodule provided by e. 

The convolution exponential exp^ : Qt — > Gt, defined by 

5*5 5* n 

exp,(5) = e + 5 + — + ■■■ + — + ■■■ , (63) 

and its inverse log^ Gt — > Qt, given by 

/ s (a-e)*(a-e) , (a-e)* n 

log* a = a-e - y - + ■■■+ -1 " - '— + • • • , 64 

2 n 

allow to identify Qt as the infinitesimal form of Gt- Convergence of these series is not an issue, because 
both 5 and a — e vanish on 1, so that only a finite number of terms survive when evaluated on a tree t. 
Because infinitesimal characters vanish on products of trees, they often only involve very simple combinatorial 
problems. Then, less trivial identities can be obtained by exponentiation. We shall illustrate this fact on 
several examples in the sequel. 

To investigate the the structure of the group of characters Gt, it is convenient to first define a decreasing 
sequence of subgroups of Gt by requiring that elements of are characters that vanish on all trees of 
order < n. This is a sequence of normal subgroups of Gt, so that one can define an increasing sequence of 
quotient groups G T = G T /K^ +1 . G T is a finite dimensional nilpotent Lie group with Lie algebra Q^ and 
one recovers Gt and its Lie algebra as the projective limits of the corresponding sequences. Besides, one 
can construct recursively these groups since G T +1 is an abelian extension of G T . Recall that given a group 
G and an abelian group K, a abelian extension of G by K is a new group G defined as follows. As a set 
G = G x K, but its product law differs from the direct product law by a 2-cocycle 9, 

(a, A) • (/?, n) = (af3, A + fx + 9(a, 0)) , (65) 

for any a,/3 e G and A, fi G K. The associativity of the multiplication on G, is ensured by the cocycle 
condition 

9{fi, 7) - 6(a0, 7) + 0(a, /3 7 ) - 9{a, (3) = (66) 

for all a, 0-f € G. 

For any rooted tree t of order n + 1, let us define 



t(a,[3)= Yl a(t 1 )---a(t n )(3[-——), (67) 



c admissible cuts of t 
c^cmpty c^full 

where we exclude the empty and the full cuts. Using the reduced coproduct A' = A — 1 ® Id — Id ® 1, the 
cocycle reads, 

9 t (o,j3) = (a0j8)oA'(i), (68) 

and thus identifies with A' considered as a multiplicative function from the algebra generated by trees of 
order n + 1 to the tensor product of to copies of the algebra generated by trees of order less or equal to n. 
This allows for a quick check of the cocycle condition since it follows directly from the associativity of A. 

TIt is graded by the number of vertices |i| and fits into the general framework of graded and commutative 
Hopf algebras. Recall that a graded and commutative Hopf algebra is a direct sum Ti — ® n( zf{H n with 
T~tm'7~tn C Ti-mn and AH n C 9meMWm®W„_ m . This framework allows us to implement inductive reasoning. 
If we further assume that Ho is one dimensional, a general structure theorem [TS] states that H. is a free 
commutative algebra so that all the previous results still hold, since most of the proofs rely on inductive 
reasoning. 

3.3 Power series of non linear operators 



When solving the differential equation (52 1 using a sum indexed by rooted trees, we associate each tree with 
a combination of the derivative of X. In order to generalize this construction, it is fruitful to view trees as 
indices for formal power series of non linear operators. 
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To proceed, let us consider any smooth map X from a Banach space £ to itself and let us define X raised 
to the power of the tree t as 



X* = Y[X^ (69) 



vet 



where n v is number of edges leaving the vertex v and X^ is the n th order differential of X . By convention, 
the tree is oriented from the root to the leaves and we compose the operators following this order. When 
evaluated at a point x € £ , X^{x) is n-linear map and the product refers to the composition of those linear 
maps from the root to the leaves. As a result, X 1 is a smooth map from £ to itself. For instance, we have 

I 

Xl = X'[X'[X]] and X<^> = X"[X,X]. (70) 
Note that this corresponds exactly to the combinations of partial derivatives used for solving the differential 



equation ( 52 ) . Let us note that there is an interesting similarity with the parenthesized words introduced 
by Kreimer in the very first paper on the Hopf algebras of renormalization |17j . 
Given an arbitrary character a £ Gt, let us define the formal power series 

* a (X) = 5>(t):J, (71) 

with St the symmetry factor of the tree, which is the cardinal of its automorphism group. This has to be 
understood as a formal power series of smooth maps, i.e. when evaluated on any x £ £ it yields a formal 
power series of elements of £, which is ordered in increasing number of vertices of the trees. Let us also 
note that the series always starts with the empty tree which is nothing but the unit of the algebra TLt- By 
convention, any operator raised to the power of the empty tree is the identity, X® = id. 

If we now consider another (3 € Gt, then we can compose ^ a (X) with ipp^X), Taylor expand and 
reorganize the series in terms of rooted trees. It turns out that the result is naturally expressed using the 
product law in Gt, 

9 a (X)o 9 fi (X) = p » a (X) . (72) 

composition convolution 

This composition law is identical to the composiion law of B-series |16j . underlying Butcher's work on the 
algebraic structure of Runge-Kutta methods. To prove this relation, we write 8 = log* (a) and consider both 
sides of as a function of s S R for a s = exp + (s6). Then, we derive both sides with respect to s and check that 
the derivatives agree using simple combinatorial identities. Let us note that if X is linear, all derivatives 
X^ vanish for n > 1 so that only trees without branchings appear. In this case eq. (72 1 amounts to the 
ordinary multiplication of power series of a linear operator. 

It is also convenient to allow the trees to carry some decorations on their vertices. Then, all the structure 
remains similar, with R, c (t) and P c (t) carrying decorations inherited from t in the definition of the coproduct. 
The main interest of decorated rooted trees lies in the fact that they allow us to define power series of several 
non linear operators, each operator being labeled by its decoration. 

The most useful series is the geometric one, defined by summing over all trees with a weight one, 

(id-X)" 1 ^^. (73) 

As suggested by the notation, it yields an inverse of id — X for the composition law. Accordingly, it provides 
us with a perturbative solution of the fixed point equation 

x = xu + X(x) (74) 

as a sum over trees 

t 



x=(id-X)-\x ) = y^(x Q ). IT:,) 



t 



S 



Moreover, let us note that the tree like structure of the geometric expansion proves to be extremely helpful 
in dealing with convergence issues. To proceed, let us assume that the n th order differential is a bounded 
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n-linear map such that ||X(")(a; )|| < a n , with a n a sequence of real numbers such that 

oo 

converges. Then, the series 



n=0 



x = (id - fcX) _1 (so) = E^'fW (77) 

is bounded in norm by the series associated to the equation z = hf(z), which is known to be an analytic 
function of h at the origin by the implicit function theorem. 

Beyond the geometric scries, it may also be useful to consider the binomial series. For any character a 
and complex number a, let us define 

/ \* n /x a(a—l), . , . a(a — 1) • • • (a — n+ 1) , 

(a)* a = e + a(a-e) + — -(a-e) * (a-e) + • • • + - L (a-e)* n + ■■■ . (78 

2 n! 

Obviously, the binomial series fulfills all the identities we expect (a)*° * (a)* b = (a)* a+b and ((a)* a )* b = 
(a)* ab at the notable exception of (a)* a *(/3)* a = (a*/3)* a , which does not hold since the convolution product 
is not commutative. 

If we take a to be the character that takes the value 1 on the tree with one vertex and vanishes on any 
other non trivial tree, then the iterations of the coproduct are easy to compute, 

(ar(t) = E Mn,0 O( °" 1) " nl ( °" W + 1) , (79) 

n=d t 

where N(n,t) is the number of surjective maps form the vertices oft to {1, • • • ,n}, strictly increasing from 
the root to the leaves, dt is the depth of the tree which is the length of the longest path from the root to 
the leaves. 

Accordingly, one can apply the binomial series to a non linear operator, 

(id + Xr = Y / J2 N(n, t) a(a - 1} --; a - Tt+1) f (80) 

t n—dt 

the coefficients in from of ^- can be considered as tree-lie generalizations of the binomial coefficients, to 
which they reduce on trees without branchings. For instance, up to order 4 we have 

(id + X) a = 1 + aX + a(a ~ 1} (x'[X] + \x"[X, X]+\ + fl(a ~ ^ (fl ~ 2) (x'[X'[X]) + X"[X, X] 

(81) 

This allows, for instance, to compute the square root of a diffeomorphism 



VidTx = -X--X'[X\ + ^X'[X'[X\]-—X'[X'[X'[X\]] 

+ ^X"[X., X'[X}} J^X'[X"[X, X}} + ^X"'[X, X,X] + --- (82) 

which fulfills \Jk\ + X o -y/id + X = id + X up to terms of fifth order in X. 

Let us also introduce the number JV(n, t) of surjective maps form the vertices of t to {1, • • • , n}, increasing 
(but not necessarily in a strict manner) from the root to the leaves. Then, by applying the binomial series 
to the geometric one, we get 

(id - at* = (da - xrr = E E %, ^ (ffi " 1) -J fl " n+1) f (83) 

t n— 1 



Then, by comparing with the expression of (id — X) a derived from (80 1, we get 

EAr(n,t) a(a ~ 1) --j a ~ n + 1) = (-1)1*1 E N(n,t)(-l) n aia + 1) "'\ a + n + 1) (84) 
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By evaluating for various values of a, we get various relations between N(n,t) and N(n,t). For example, if 
a = — 1 we obtain the sum rule 

1*1 

i=(-i)i ( i^(-ir%t), (85) 

n=d t 

which is a convenient way to check the computation of N(n, i). For a = —2, we get 

1*1 

2 + N(2,t) = £(-l)"+l'l(n + (86) 

n=d t 

which counts the number of terms in the coproduct. 

The geometric series and its cousins can be used in two complementary ways. In many cases, the 



equations we encounter in physics can be written as fixed point equations of the type (74 1, so that rooted 
trees allow for an efficient determination of their perturbative solutions. For instance, this is the case for 
the self-consistency condition derived in mean field theory or the Schwiger-Dyson equations in QFT. In 
mathematics, this applies to the solution of an ordinary differential equation written in integral form, as we 



shall work out in section 3.5 or to the Lagrange inversion formula. 



Alternatively, one can be given a pcrtubativc expansion exhibiting a natural tree like structure. If one can 



find an operator X such that this expansion can be written as a geometric series ( 73 1 , then this expansion is 



known to be a solution of the fixed point equation ( 74 ) . Accordingly, the solution of this equation provides 



a good candidate for the sum of the tree like expansion. We illustrate this procedure in the next section. 
3.4 Wigner's semi-circle law from rooted trees 

Let us give a very simple illustration of way one can use rooted trees to recover Wigner's semi-circle law of an 
ensemble of GauBian Hermitian N x N matrices. This technique is nothing but a rooted tree formalisation 
of the treatment that can be found in [6] . 

We are interested in the large N behavior of the density of eigenvalues, defined as 

„(A) = Im^ J [DM] e -™ ( M) £ S{X_X^ (g?) 

i 

with Xi denoting the eigenvalues of M and V(M) — \M 2 . [DM] — \\ i dMu Yli<j dReMijdlmMij is standard 
Lebesgue measure on the space of N x N Hermitian matrices and S the Dirac distribution. 
Using 6(x) = ^Im (^^zj^j with e — ► + , we can rewrite p as 

Introducing z = A — ie to simplify the notation, we expand the resolvent as a geometrical series 

' - + -M- + -M-M- + ■■■ . (89) 



z — M z z z z z z 

The expectation value G(z) of z ^ M is easily computed using Wick's theorem with the GauBian matrix 
integral. It leads to Feynman rules with a single line propagator for ^ a double line propagator for the 
propagation of the matrix field and a vertex set to 1 for the interaction of 2 single lines with a double line. 
Then, each face of the resulting diagram yields a factor N because of the summation over matrix indices. 
Thus, the leading contribution in the large N limit is given by those diagrams that have as many faces as 
possible for a fixed number of propagator. These diagrams are nothing but the celebrated planar diagrams. 

The planar Feynman diagrams that contribute to G(z) are in one to one correspondence with planar 
rooted trees (i.e. rooted trees counted as different if they differ by a permutation of their vertices). The 
bijection is obtained by drawing boxes that are either disjoint or nested around the matrix propagators, with 
the diagram always in the outermost box, even if there is no matrix propagator. This outermost box defines 
the root and the tree is drawn from the root to the leaves by following the natural hierarchal structure of 
the inclusion of the boxes. 
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Figure 7: A planar diagram contributing to G{z) 



Isomorphic planar trees correspond to Feynman diagram having the same contribution to G{z) . Moreover, 
recall that a rooted tree is an isomorphism class of planar rooted tree and the sum over all representatives 

with 

1 



X 1 ( 0) 

of the class of t equals g A, ; , 



X{x) 



(90) 



Note that the n th order derivative of X generates a factor nl at each vertex with n outgoing edges. When 
multiplied together along the tree and divided by the symmetry factor, the resulting number counts the 
cardinality of the isomorphism class of planar trees corresponding to t. 
Therefore, G(z) is the sum of the non linear geometric series 



rM _ V It(0) T 



(91) 



with In the N x N identity matrix. The sum of the geometric series is the solution x of the quadratic fixed 
point equation X(x) — x with leading term ^, so that 



G{z) 



In- 



(92) 



and we recover Wigner's semicircle law as 



p(A) = ^V4^, 



with support in [—2, +2]. 

If we expand G(z) in powers of -, we recover the ubiquitous Catalan numbers C, — — 



(n!) 2 (n+l) ' 



oo 



n=0 



(93) 



(94) 



This is not surprising since the Catalan number C n count, amongst many other things |18j . the number of non 
crossing partitions of a set of 2n elements into blocks of 2. More sophisticated examples of summation over 
rooted trees applied to random matrices are described in the lecture by P. Di Francesco 19J. In particular, 
it is shown how some similar techniques can be used for a quartic potential V(M) 



|M 4 . 



3.5 Perturbative solution of differential equation 

To solve pertubatively a time dependent non linear differential equation 



dx 
ds 



X s (x), 



(95) 
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Figure 8: Geometric interpretation of the coproduct 



with boundary condition x(s ) — x , it is convenient to write it in integral form, 

x(s)=x + [ ds' X s ,(x(s')). (96) 



This is a fixed point equation for the integral non linear operator x(s) i— > ds' X(x(s')) with the function 
s i — ► x(s) as the unknown. Using the geometric series, it solution is expanded over trees as 

^ s )=£^ L d W sf[x^\x ). (97) 



vet 



To explain the meaning of these terms, recall that any tree is orientated form the root to the leaves. Then, 
associate to each vertex v a real variable s v in such a way that sq < s v < s v+ if v + is a vertex immediately 
above v and s v+ = s for the root. This defines what we call a "treeplex" P s SQ C K'*', that reduces to a 

ordinary simplex if t is a tree without branching. Then compose all the multilinear operators X^ v \xo) from 
the root to the leaves and integrate over the variables s v . Obviously, for linear X we recover the ordinary 
time ordered products, with an integral over a simplex instead of a treeplex. 



If we write the solution of (95 1, as f2 SiSo (xo) = ^ t £l s so ( x o)i ^ ne semigroup property holds 



admissible cuts c 



where to is the part of t that contains its root after the cut and t\, . . . ,t n are the branches we cut. The 
symbol o means composition followed by Taylor expansion and o c is the same operation, to be performed 
at each vertex immediately above the cut. This relation is similar to the definition of the coproduct and an 
analogous equation exists for the antipode, related to the inversion of fi s ,s ■ 



Comparing both sides of (98 1 yields a geometric interpretation of the coproduct: Cut P by a plane 



Sroot — s ii so that it is written as a disjoint union 

J Lso= U 6c(^, So x---x^ o x^ Si ) (99) 

C admissible cut 

with P c (t) = ii . . .t n and R c (t) = to. 6 C a suitable permutation of the labels preserving the ordering of the 
tree. For example, for the simplest tree with one branching, this is illustrated on figure [8] and the polytopes 



in this decomposition are in direct correspondence with the terms in ( 58 ) 



For a time independent equation, the integral over 7* So factorizes and defines the tree factorial 



d |4| s=i. (100) 
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The latter can be computed recursively: If t is a tree obtained by grafting the trees t\,...,t n to the root, 
then 

f! = (|ti| + "- + M + l) til- ••*»!• (101) 
For a tree without branching we recover the ordinary factorial t\ = \t\\. 

Accordingly, the map Q s ,s is a power series of the non linear operator X 

1 (s-so) 1 * 1 



(1 ° 2) 

which is a non linear analogue of the exponential. The exponential character can be written as 

exp, {(* - s )S} (t) = (5 ~' o)|t ' (103) 

with d € Qt the infinitesimal character that takes the value 1 on the simplest non empty tree and vanishes 
otherwise. It satisfies the linear equation 

^- exp,, {(s - s )8} = 5* exp„ {(s - s )8} = exp* {(a - s )^} * 6. (104) 
as 

Thus, we have traded a non linear equation for a linear one, at the price of working with power series of non 
linear operators. 

The relation ( 103 I provides a combinatorial definition of the tree factorial, in analogy with what we did 
for the binomial series (78 1. Because 5 is an infinitesimal character, there arc tremendous simplifications in 
the convolution powers of 6 are evaluated on the iterated coproduct and only one term survives 

6* n (t) = N(n,t) (105) 

if n = \t\ and vanishes otherwise. Therefore, ^ — N ^'^ so that 

^ # {labelings of the set of vertices in {1, ... , n}} (106) 

# {labelings of the set of vertices decreasing from the root to the leaves} 



Turning back to numerical analysis, the Runge-Kutta method (55) one naturally associates a character 
a constructed as follows. Decorate the vertices of the tree t by indices in {1,2, . . . ,n} and associate each 
edge from i to j with the matrix clement a^ . Multiply all these matrix elements together with a factor bi 
associated to the root and sum over all indices. Then, the approximate solution after one step reads 

x 1 = J2~a(t)X t (x ). (107) 

It is a method of order n if a coincides with the exponential on all trees of order at most n. This defines 
another formal power series of non linear operators and the composition law of Runge-Kutta methods is 
represented by the convolution of their characters. 



4 Renormalization, effective actions and Feynman diagrams 
4.1 Renormalization for a differential equation 

Before dealing with path integrals and Feynman diagrams, it is fruitful to first investigate the perturbative 
renormalization for an ordinary differential equation [20 . Recall that the differential renormalization group 



equation (46) reads 

A^ = /9(A,S A ), (108) 

with boundary condition S\ — Sq at a very high energy scale Ao. In physics, 5a is an effective action that 
encodes the description of the physical system under scrutiny, valid up to an energy scale A, while So is 
the bare action which agrees with the effective action at the cut-off scale Ao- Ignoring temporarily the path 
integral interpretation of the renormalization group equation, let us simply consider a general equation of 



the form ( 108 1 for a variable 5 in a Banach space £. Following the general interpretation of renormalization, 
we are interested in the low energy behavior of 5a for A <C Ao. 

Because all the physical quantities can be measured in unit of massj^] we further assume that the space £ 



3 In the system of units h = c = 1 which we adopt here. 
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comes equipped with an action of the dilatation group This action is written in exponential form as e sT> , 
with T> a linear operator. It is often convenient to assume that T> is diagonal, with eigenvalues corresponding 
to the canonical dimensions. Dimensional analysis is assumed to be compatible with the renormalisation 
group equation, 

/3(e s A, e sV S) = e sV [3(A, S). (109) 

To take dimensional analysis into account, it is convenient to introduce the dimcnsionlcss variables s = 
log(A/A ref ) with A rof a fixed reference scale, u(s) = A~ v -S\ and X(u) = A~jP -/3(A rof , <Sa). In terms of these 
new variables, the renormalization group differential equation takes the simpler form 

fill 

^ = -Vu + X(u), (110) 

with boundary condition u(s ) = u . 

The pertubative solution of this equation is easily obtained by writing it in integral form 

u(a) = e-^-^uo + f dtfe<'-'"> v X{u{s')). (Ill) 

Making use of the geometric series of the integral operator T ■ u(s) = J° ds' c( s ~ s ~> v X(u(s')), 

u = (id - T)- 1 (u ) (112) 

we expand u(s) as a sum over trees, using the following rules to compute the contribution of each tree t. 
First, orient all the edges from the root to the leaves and then associate 

• variables s v to the vertices that leave in the the treeplex J* s . 

• linear operators e~( s "+~ Sv -' > ' D to edges from v + to v_ (recall that edges are orented from the root to 
the leaves); 

• multilinear operators X( nv \e~( Sv ~ So ^ v Uo) to vertices with n v outgoing edges. 

Finally, compose from the root to the leaves all operators associated to edges and vertices, multiply on 
the left by c~( s ~ Sroot )' 5 and integrate over the variables associated to the vertices. This is very similar to the 
expansion we encountered in the previous section for a time dependent equation, at the notable exception 
of the insertion of the of the linear operators e^ s "+ ~ s "- ^ v along the edges. 

To simplify the notations, let us abbreviate as X$ the n-linear map obtained by differentiating X n 
times at e~^ Sqj ^ So ^ v u . For example, the linear tree of order two contributes as 

fds 2 [ S2 d Sl e-('-»Wx'Je-(»- ai '> 7 >X ai ], (113) 

whereas the order 3 tree with one branching yields 

\ f "ds 3 r ds 2 r d Sl e-^ v X^e-^-^ v X S2 ,c-^-^ v X Sl ] (114) 

^ J So J s J s 

In the simplest case of a linear X, the contributions of all branched trees vanish and we recover Duhamel's 
formula (or in the physics terminology, the Schwingcr-dc Witt expansion) of the exponential of —T> + X in 
powers of X, 

In analogy with pcrturbativc renormalization in QFT, let us try to get rid of the cut-off by taking the 
limit so ~ > 00 • If we discard all the non linear terms, the analysis us easy to perform: 

• variables u 1 " that belong to the negative eigenspaces of V (irrelevant variables) tend to 0; 

• variables w mar that belong to the zero eigenstpace of V (marginal variables) remain constant; 

• variables w rcl that belong to the positive eigenspaces of V (relevant variables) grow to infinity. 



23 



In perturbation theory, this picture is slightly modified by the addition of the non linear terms. Marginal 
variables have a polynomial dependence on s, analogous to the log A terms in QFT. Moreover, irrelevant 
terms at a given order depend on the marginal and relevant ones at a lower order, so that they also experience 
divergences. This phenomena is analogous to the subdivergences of QFT. 

The solution of the problem of UV divergencies is natural from the Wilsonian viewpoint: Instead of 
imposing the boundary condition for the divergent variables at the very high energy sq tending to infinity, 
impose the latter at a fixed low energy s r . This corresponds to a determination of the relevant (and marginal) 
parameters by a measurement at low energy. To implement this change of boundary condition, let us denote 
by V the projector onto the relevant and marginal variables and by T' the non linear operator 

T'[u\{s) = V I'" da' e (s - s ' )v X(u(s')). (115) 



Therefore, the renormalized theory, expressed as a function of the renormalized parameters (irrelevant para- 
maters at so, relevant and marginal ones at s r ), is obtained as a fixed point for the renormalized operator 
T + T' . To compare the renormalized and unrenormalized expansions, let us write 

1 — 1 r 1-1 r , r l 

\&-{T + T') = SA.-T o id - T' o (id - T) . (116) 

In the QFT language, the LHS is the renormalized theory while the RHS is the unrenormalized one composed 
with the counterterms. The latter are the boundary conditions on the relevant and marginal parameters 
to be imposed at the high energy sq in such a way that these parameters match their fixed values at low 
energy s r . At the level of the tree expansions, the expansion of the renormalized theory can be pictured 
as a sum over trees with both white vertices for T and black ones for J 7 ', while the unrenormalized theory 
only involves white vertices. Counterterms only involve trees with a black root, followed by black and white 
vertices that correspond to the renormalization of the subdivergences. This is the analogue of the BPHZ 
formula for a differential equation. 

In the renormalized expansion, all the so dependence has dropped from the divergent exponentials and 
their integrals, so that the contribution of each tree admits a finite limit as s — > oo.. Moreover, the initial 
conditions u n ' for the irrelevant paremeters all enter through the combination e~ x '( s ~ ;i0 )it rr , which goes to 
as so — * oo- I n this way we recover universality: The renormalized theory does not dependent on the choice 
of the irrelevant parameters at high energy. 

When formulated in terms of differential equations, the renormalization procedure is general enough to 



encompass Polchinski's equation (49 1, as far as perturbation is concerned. Solving the latter using rooted 
trees yields a smart of proof of the rcnormalizability of theory in D = 4 that does not refer to Feynman 
diagrams [ST] , in addition to the original proof by Polchinski [TO] . 

Though rooted trees are a convenient device that often leads to convergent expansions, it is fair to say that 
this is seldom the case for differential equation with boundary condition at oo. To illustrate the issue raised 
by non perturbative computations within the differential renormalization group approach, let us consider a 
very simple example with one marginal and one irrelevant coupling, 

■ * = /V (v = o), 

I = u + 1 g (2?=-l). 1 n 

This system can be solved using the pertubative techniques we have presented. However, because of the 
extreme simplicity of these equations, it is much more convenient to determine an exact solution. Let us 
first investigate the case of the marginal coupling constant g. In the unrenormalized case, the boundary 
condition is imposed at very high energy, so that the unrenormalized solution and its expansion read 

oo 

9(s) = r = .90 J2^o) n (s - s )". (118) 

1 - f3g (s - s ) ^ 

In terms of rooted trees, the LHS is the sum over rooted trees with vertices having only at most two outgoing 
edges (since the operator g i— > (3g 2 is quadratic) and a weight corresponding to the tree factorial. 

As expected, the coefficients of the expansion, which play the role of the Feynman diagrams, diverge as 
so — > oo at fixed go- Perturbative renormalization is achieved by imposing boundary conditions at a fixed 
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g 



boundary condition g(s r )-g r 




s = 15 



Figure 9: Convergence towards the renormalized theory for the differential equation (1171 



low energy scale <?(s r ) = g r , so that the rcnormalization coupling constant reads 



9(s) = 



1 - f3g T (s - s r ) 



(119) 



n=0 



The transition from the unrenormalized coupling to renormalized one is achieved by choosing go in such a 
way that g obeys the new boundary condition g(s r ) — <? r , which is easily seen to be 



90 



i - PgAso - Sr) 



(120) 



n=0 



Then, as we substitute this expresssion of go in (118) we recover (1191. At the perturbative level, this 
substitution is best seen using rooted trees and illustrates the way BPHZ rcnormalization must be performed 
with the trees. 



Once the renormalized coupling constant ( 119 1 is known, the solution of the second equation follows from 
the method of variation of constants, 



-(s'-s) 



u(s) = u Q e s - s °+ ds'- 

J so 1 - P9x{s' ~ Sr) 



(121) 



First, let us notice that the u is independent of uo as so — > oo. This is universality in this simple context: 
The renormalized theory does not depend on the irrelevant terms at high energy. At s r — and so — ► oo, 
the exact solution can be expanded as 



u(9t) = -7 / dse s 9l 

Jo 1 - 9vPS 



J20 n nl(gr) n+1 . 



(122) 



n=0 



Though it provides accurate approximations of u for small g ri it is a divergent power series, Borel summable 
only if (3 < 0. Nevertheless, if /3 < the coupling g(s) remains bounded and u is well defined as a functional 
of g : g(s), 



r 

H = / 

Jo 



dse s g(s) 



(123) 



This is a particularly simple example of an effective expansion |llj . whose counterpart in QFT is the fact 
that a asymptotically free theory, characterized by (3 < 0, may be well defined when expanded in a scale 



dependent coupling instead g(s) of a single renormalized coupling g r . In case (3 > 0, (121 1 can still be given 
a mathematical meaning by deforming the contour of integration in the s plane. However, in this case there 
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axe extra ambiguities related to the monodromy around the pole so that the construction of the theory 
requires extra information not contained in the perturbation theory. Nevertheless, the theory may lose its 
physical meaning because it is no longer real. 

The occurrence of the previous type of n! behavior that prevents a naive summation of the perturbative 
solution with boundary conditions at infinity is due to the phenomenon of resonances. Consider a general 
renormalization group equation 

du „ , . , 

— = -Vu + X{u) (124) 

and assume that T> has been diagonalized in a basis (e^) with eigenvalues A^. We further assume that X can 
be written as a polynomial (or an analytic function) 



x( U )= ^i.. jn (u h r^---(u jn r-. (125) 



A monomial in (125) is resonant if there is a non trivial relation of the type 

Xi — m 1 \j 1 + ■ ■ ■ + m n Xj n , (126) 

with rrik £ N. 

As we try to solve the renormalization group equation pertubatively in u, the Poincare-Dulac theorem 
states that, by a series of successive changes of variables of the type v = u + F(u) with F having at least a 
second order term, the renormalization group equation takes the simpler form 

d ^- = -Vv + Y{v) (127) 



with only resonant monomials in Y . Because the resonance condition ( 126 ) can be formulated as 

e~ sI> oY = Y o e~ sV , (128) 



all the exponentials but the the first one cancel in the perturbative expansion of the solution given in (112) 
Therefore, the iterated integral along the tree generates a power law s n and leads to a factorial behavior n 
when integrated together with the exponential along the root. 

Let us conclude this section by two remarks that may be the subject of some future work. 



• It could be interesting to derive the normal form of the QFT renormalization group equation ( 49 ) , 
with a view towards a non perturbative analysis. This could be done along the lines initiated in [22 
for a toy model. 

• Reduction to normal forms is a non linear generalization of the Lie algebra decomposition into com- 
muting semi-simple and nilpotent elements |23j . As such, it could be conveniently described in the 
algebraic framework presented in the sections devoted to rooted trees and Feynman diagrams. 



4.2 Hopf algebra of Feynman diagrams and effective actions 

In the view of the algebraic structure we have presented in the previous sections, let us investigate Wilsonian 
renormalization in QFT. The basic renormalization group transform is the integration over fast modes that 
induces the transformation of effective actions S — > S', 

e s'm = tf J [D X ] e -sx- A_1 -x e 5 ^+*l. (129) 

A is a suitable propagator that implements the integration over momenta lying between A and Ao. In what 
followos, we consider A to be arbitrary propagator, with suitable UV and IR cut-offs. JV is a normalization 
factor, independent of S and <fi, that ensures that S' = for 5 = 0. In the sequel, we do not display such 
normalization factors anymore. Note that we have used S instead of —S in the exponential to avoid some 
irrelevant signs. 
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The rcnormalization group transformation of effective actions is best formulated using the background 
field method. The rationale of the background field method is to split the field into a fixed classical field <t> 
and a quantum field X over which we integrate, so that we end up with an effective action function for the 
classical field (p. In our setting, the fast modes play the role of the quantum field whereas the low energy 
ones remain as an argument of the low energy effective action. Using this technique, S' is expanded over 
connected vacuum Feynman diagrams (i.e. diagrams without external leg) as 

S'M = £^%]. (130) 

7 7 

The sum starts with the trivial empty diagram whose contribution is S [4>\. For other diagrams, A 1 {S) 
is computed using the background field technique Feynman rules: Propagators are determined by A and 
n-valent vertices are n th order functional differentials of S at cf). As is always the case for diagrammatic 
expansions, one has to further divide by the symmetry factor S 7 , which is nothing but the cardinal of the 
automorphism group of the diagram. All the dependence on <f> is encoded in the vertices so that there is no 
need for external sources and we restrict ourselves to vacuum diagrams. Note that there is no restriction 
on S since the convergence of each Feynman diagram is ensured by the cut-off. Strictly speaking, this true 
only for a sharp cut-off, when considering a cut-off with a exponential decrease we have to assume that the 
momentum dependence of the functional derivatives of S does not grow faster than a polynomial. 

To illustrate the background field technique, let us give a simple example with a finite dimensional integral 
and a positive definite symmetric matrix Aij as a propagator. 

X^^X ±_ ST d 3 S d 3 S 

"■ - .. 12.4^. d X ^d X ^d X ^ W <lJl l2 ' n i3 ' h d X ^d X ^X h[ ) 



In the spirit of Wilsonian renormalization, one has to pursue this procedure and consider S' to be a 
starting point for a new integration over yet slower modes. This is the reason why we are forced to deal 
with actions of a very general form. Indeed, even if S is a simple </> 4 interaction, this is no longer the case 
for S' that contains arbitrarily complicated vertices. To handle this type of computation, it is necessary 
to modify the Hopf algebra defined by Connes and Kreimer as follows. Let Ti-p be the free commutative 
algebra generated by all connected vacuum Feynman diagrams with vertices of arbitrary valence, including 
univalent and bivalent ones. It is equipped with a Hopf algebra structure with trivial counit and coproduct 

A( 7 ) = 7 ®1 + 1®7+ 7i"-7n®7/(7i"-7n), (131) 

71>--->7n 
disjoint subdiagrams 

where the sum runs over all possible mutually disjoint connected subdiagrams. In our setting, a subdiagram 
is a subset lines and vertices related to them. The reduced diagram r/( 7 i • • -j n ) is obtained by shrinking 
each 7^ to a single vertex in T. Obviously, Hp is graded by the number of internal lines I but not by 
the number of vertices. However, the number of loops L also provides a grading because of the identity 
L = I — F + 1 for a connected diagram, but this grading does not yield a connected algebra since all the 
trees have grading 0. To these gradings is associated the two parameter group of automorphisms 

<Pa,bCr) = a' 7 ^ 7 7, (132) 

with Ij (resp. L~) being the number of internal lines (resp. loops) of 7. Because of the grading, one can 
apply the general theory of graded and commutative Hopf algebras: The antipode may be defined recursively 
and characters of Tip form a group Gf, with a Lie algebra Qp of infinitesimal characters. Note that elements 
of both Gf and Qp are uniquely defined by their values on the diagrams, which may assume any values in 
a given commutative algebra. As for trees, the exponential defines a map exp,, : Qp — > Gf by 

6*5 6* n 

exp,(5) = e + 5+ — + ... + _ + ... , (133) 
which is a solution of the differential equation 

^ cx P>t (s<5) = 6 * exp t (s5) = exp^(sS) * 6 (134) 
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reducing to the unit of Gt at s = 0. 

The elements of Gf can be considered as weights of the for the Feynman diagrams in the expansion of 
the transformation of effective actions : S — > S' 

S' = $^(S)=2«(7)^ (135) 

7 ' 

Given any two characters a, (3 £ Gf, there is a composition law analogous to the one we encounter for power 



*fl"^= <*a ■ (136) 



series of non linear operators in (|72|, 

composition convolution 

This composition law describes how the weight of the Feynman diagrams entering in an iteration of two 
effective action computations can be obtained by by inserting diagrams with weight a into diagrams with 
weight (3. To prove this relation, we introduce 5 = log* (a) and a s = exp < ,(s5). At first order in s, this relation 
does not involve multiple insertions of diagrams and is easily verified by simple combinatorics. Then, we 
integrate the differential equation to recover ( p6| . In a certain sense, it is a Wilsonian analogue of the 



morphism from the group of graphs into the diffeomorphisms of the coupling constant presented in [3] . 

In physics, there are many characters that appear naturally. First, one may consider a to be the character 
that selects 1PI diagrams and (3 the one the selects the trees. Then, a * (3 takes the value one on all 



diagrams and the composition law ( 136 1 simply states that the one recovers the complete effective action 
by first computing 1PI diagrams and then inserting into trees. As another simple example, let us consider 
the character a s ,h{l) — s I ~ih l ~< , so that the corresponding effective action transformation S — + S' can be 
expressed as 

S'[<f>] = hlog{ /[£>*] e-**-< fcsA > _1 -*e^V (137) 



As we expand 5" in term of Feynman diagrams, we find a weight hs for every internal line and hr x for every 
vertex together with a global factor of h. Thus a connected diagram with 1 1 internal lines, V 7 vertices and 
L.y = Ij — V^ + 1 loops is weighted by s Il h I ~<~ v ~> +1 — s I ~'h L ~< . Other characters of interest are the polynomials 
in N that arise in the O(N) invariant models. 

To handle more complicated situations, it may necessary to decorate the Feynman diagrams. For instance, 
if we want to perform the second effective action with a different propagator B, then we have to decorate 
the internal lines with a index that differentiates A and B. With decorations on the internal lines, the Hopf 



algebra structure remains similar with the subgraph that we extract in ( 131 1 carrying the decorations. Then, 
we can perform a first effective action computation with propagator A and character a, followed by a second 
one with B and (3. The results amounts to a single computation with decorated internal lines that selects 



the type of propagators we use. Then, the analogue of ( 136 ) reads 



*f°^=C (138) 

Both a and [3 are characters of the Hopf algebra with two decorations on the internal lines, but a and (resp. 
(3) vanish on diagrams that do not contain A (resp. B). The previous relation is a generalization of the 
convolution of Gaufiian integrals 

j[D X ] j[Dt]vv{-\xA- l x- ldB- 1 Z + S[ct> + X + Z}} = 

J [D(] exp {-|C(A + 5)"^ + + C]}, (139) 



which we obtain from ( 138 1 when the characters a and (3 only take the values or 1, depending on whether 
the diagrams only involve the corresponding propagator. Obviously, this construction can be generalized to 
an arbitrary number of propagators so that we can iterate the effective action construction. Unfortunately, 
this construction is not sufficient to handle Polchinski's equation since the latter involve an infinite number 
of iterations corresponding to the differential equation. Moreover, the various types interactions (relevant, 
marginal and irrelevant) play fundamentally different roles so that it is necessary to distinguish them in the 
Feynman expansion. We postpone these questions to the section |4.4| devoted to the iterative solution of 
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Polchinski's equation. Before we deal with the latter, let us investigate some algebraic properties of the TLp 
and illustrate them in the next section on some properties of the Tutte polynomial. 

First of all, let us notice that there are only two diagrams of the degree one (i.e. with one internal 
line): The bridge with two vertices and the self- loop with a single vertex. Let us define two inifinitesimal 
characters: 5 trcc takes the value one on the bridge and vanishes otherwise, and (hoop takes the value one on 
the self- loop and vanishes otherwise, 

= 1 -d f WO) = J m 



KUl) = if I Ko V {l) = if 7 ^0 

By iterating the coproduct, it is easy to see that 

exp, {aS tICC + b 6i oop } (7) = a I ^ L -<b L \ (141) 

where 7 7 is the number of internal lines and L~ t the number of loops. By coassociativity, the iteration can 
be performed either on the left or on the right. If we perform it on the left, then we contract the internal 
lines and we get a factor of a if the line is a self-loop or b if it is not. When performed on the right we delete 
internal lines and we also get factors of a and b depending on whether the line is a bridge (i.e. a line whose 
removal disconnects the diagram) or not. Note that all contraction/deletion schemes yield the same result 
and can be mixed together: At one step we can delete and then contract at another. 

Any infinitesimal character 8 G Qf defines two derivations of Hf that correspond to the infinitesimal 
form of right and left multiplication by exp^(sS) on functions over Gt at first order in s , 

f<8 = (<5®Id)oA(/), 

8\>f = (Id<8>£)oA(/), 1 > 

for any / € Ht- Obviously, left and right action for different characters commute. In case of the infinitesimal 
characters 8 tlee and 8\ oop , these operations have natural interpretations in terms of deletions and contraction 
of edges. Recall that a self-loop is an edge with both ends linked to the same vertex and a bridge is an 
edge whose removal disconnects the diagram. Note that an edge connecting a single vertex to the rest of the 
diagram is considered as a bridge. 

• <5i oop E> 7 is a sum over all the diagrams obtained from 7 by cutting a line which is not a bridge. 

• <5tree > 7 is a sum over all the diagrams obtained from 7 by cutting a bridge. 

• 7 < 5ioop is a sum over all the diagrams obtained from 7 by contracting a self-loop with one edge. 



7 < 8 trC e is a sum over all the diagrams obtained from 7 by contracting a line which is not a self-loop. 



As will be illustrated in the next section, these rules allow to translate simple diagrammatic manipulations 
into differential equations. 



4.3 An interlude on the Tutte polynomial 

The Tutte polynomial is a two variable polynomial associated to graphs that plays a fundamental role in 
graph theory and statistical mechanics (see |H] for a recent review). In what follows, we shall illustrate 
its universality property in graph theory and its relation to the Potts model partition function using the 
techniques we introduced in the previous section. 

The Tutte polynomial is a polynomial in x and y which is multiplicative on disjoint unions and that can 
be defined using the rank-nullity formula on connected graphs, 

P TMe (x,y, 1 )= J2(y-ir (A) (x-l) r(E) - r(A) , (143) 
ace 

where the sum runs over all subsets of the set of edges E of 7. In the QFT language, the nullity and the 
rank of a connected diagram can be expressed in terms of the number of internal lines and loops 

fj = ft T (144) 



29 



Moreover, subsets of the sets of edges are in one-to-one correspondence with possibly disconnected subdia- 
grams, so that the Tutte polynomial can be written as the evaluation at s = 1 of the character 

a = exp,, s {Sttee + (y - 1) <Sioo P } * exp,, s {(x - 1) S tICC + <5i oop } , (145) 



where we have used ( 144 1 and ( 141 1 to write the rank and nullity using the exponentials. From this expression, 



it is obvious that the Tutte polynomial character can be expressed as 

Tutte polynomial character at x=0 on the subgraphs 
, ~ ^ 

a = exp^ s {Atcc + (y - 1) <5ioo P } * exp,, s {-8 tICC + 6\ oop } 
exp„ s{S t rcc ^loop 

} * exp, S {(x - 1) 6 t ree + ^loop}, 

v v ' 

Tutte polynomial character at y— on the reduced graphs 

which is nothing but the convolution identity of Kook, Reiner and Stanton (see [25). 



(146) 



Then, using (134 1, a can be characterized as the solution of the differential equation 
dct 

— = xa * Stree + D ^loop *a + [<5 tree ,a] t - [S\ oop , a]„ (147) 

that reduces to the trivial character at s = 0. 

This characterisation of the Tutte polynomial as the solution of a differential equation can be used 
to derive its universality property. The latter is formulated as follows. Consider a four variables graph 
polynomial Q(x,y,a,b,"/) that is multiplicative on disjoint unions and one vertex unions (i.e. diagrams 
joined at a single vertex), and has the following behavior under contraction and deletion of an arbitrary 
internal line: 

• if 7' is obtained by deleting in 7 a bridge, then 

Q(x,y,a,b,j) = x Q(x,y,a,b,j'), (148) 

• if 7' is obtained by contracting in 7 a self-loop, then 

Q(x,y,a,b,i) = yQ(x,y,a,b,i), (149) 

• if 7' is obtained by deleting in 7 a regular line (i.e. a line that is neither a bridge nor a self-loop) and 
7" obtained by contracting the same line, then 

Q(x, y, a, 6,7) = aQ(x,y,a,b,j') + bQ(x,y,a,b,j"). (150) 
Then, the character defined by 

0(7) = a I TQ{x,y,a,b,i) (151) 

obeys the differential equation 

d(3 

— = x/3 * 5 tiee + yS loop */3 + a [S tiee ,0\ t - b [S loopi [3]^ . (152) 

This follows from the discussion at the end of last section, where we have formulated algebraically the 
contraction/deletion of internal lines. Note that the contraction of a self-loop cancels in the fourth term 
while the deletion of a bridge cancel in the third term because of the multiplicativity on one vertex unions. 



Acting with the automorphism </? a -if,-i defined in ( 132 ), we obtain the Tutte polynomial differential equation 
( 147 1 with modified parameters ^ and | 



d x y 

-T^Pa-\b-i{l3) = ~ <Pa-\b-i(/3) * <5trcc + T Sloop * lfi a -i b -i(/3) 

ds a b 

+ [6t Tee ,<p a -i b -i(/3)] m ~ [<5loop, ^a-ifc" 1 (/?)]* • ( 153 ) 

Accordingly, (3 = (p a j,{a), with a the corresponding Tutte polynomial character, so that 

Q(x, y, a, 6, 7) = a^V^P^J-, f , 7). (154) 

a 
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This is the universality property of the Tutte polynomial: Any graph polynomial multiplicative on disjoint 
and one vertex unions satisfying some simple deletion/contraction identities can be expressed using the Tutte 
polynomial. 

As a second illustration, let us derive the relation of the Tutte polynomial with the g-state Potts model 
partition function using Feynman diagrams and effective actions. To this aim, we write the Tutte polynomial 
as a product of two characters ( 145 1 that we interpret as a composition of two background field effective 
action computation using ( 136 1. To simplify the notations, we set u = x — 1, v = y — 1 and q = uv, so that 
the factorization of the Tutte polynomial character reads, at s = 1, 

a = exp^ {S tICC + v <S loop } * exp,, {u S tree + S loop } . (155) 

Besides, for our purpose it is only necessary to consider effective actions for a single dimensional field which 
is nothing but a real variable. The first character is exp^ {<5t roc + v ^loop} and its induces a transformation 
of effective actions S S' that weights the Feynman diagrams with the character a(j) — v L ~< . Then, using 
( 137 1 we get 

S'[cf>} = vlog { f[D X ]e-^ 2 e^l. (156) 



Next, we consider S' as a starting point for a new effective action computation S' — > S" with a weight 
exp + {u5 tICC + <5ioop}i which amounts to weight internal lines by u and loops by ~, so that 



* log I J [D<f>] e~^ 2 e uSl ^ ] | . (157) 



As we replace S' by its expression in terms of S in the last equation, we get 



According to the general rule of effective action composition (1361, this amounts to a single effective action 
computation with Feynman diagrams weighted by their Tutte polynomial. In particular, if we take a universal 
action like S[4>] — Ae^ and evaluate S" at ip = 0, we get the Tutte polynomial multiplied by a power of A 

l -\ogU[D4>]e-^ 2 (J[D X ]e-^\^ +x \\= ]T ^P™„(u + 1, v + 1, 7 ), (159) 

7 connected ^ 

with Vy the number of vertices of 7. The choice of the universal action S[<f>] — Ae"^ is motivated by the fact 
that it generates all diagrams with a factor 1. If restrict ourselves to a monomial like A ^7-, then we generate 
only the diagrams with n-valent vertices. 

Unfortunately, it is not yet possible to interpret ( 159 1 as a generating function for the Tutte polynomial, 
since this would require us to be able to identify the contribution of each diagram in the perturbative 
expansion. To proceed, we have to write the integrals over <fi and \ 85 a single integral which is a perturbation 
of a Gaufiian integral. This can be easily done when q is an integer by introducing q independent fields x 

[ [ Dx ] e -^x 2 e is[0 +x ] V = [ H [D Xi ] e-^^^ 2 e^ s ^+^ (160) 

^ 'J l<i<q 



Then, the integral over <j> and \ m ( 159 1 reads 



II Me-^ 2 e-^^ (Xi)2 e«^ S[Xi+0] (161) 



l<i<q 

It is convenient to trade \i f° r d — Xi + 4> so that the integral over <p is Gaufiian 



J[Mj II [Dtte-KWf+tett* xe-^^« xe^iW. (162) 



l<i<q 
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Performing the GauBian integral over <p and discarding an overall constant which is absorbed in the normal- 
ization, we get 

/ ]Jm}e-M(^ 2 -^(^) 2 )} Xe i^< (163) 



l<i<q 



This can be written as a perturbed GauBian integral over a multiplct of fields £ = (£, ), so that (1591 reads 
ilogjy^Je-^^e^)^ ]T ^P Tutto ( W +l, W +l, 7 ), (164) 

7 connected ^ 

with a q x q propagator 

A = v + M (165) 
where M is the q x q matrix whose entries are all equal to 1 and an interaction 

V(0= A $> X P^)- (166) 



The expansion of (164) over Feynman diagrams can be considered as a generating functions for the Tuttc 
polynomials, since each diagram is weighted by -P Tu tto(u + lj ^ + 1,7). Moreover, the Feynman rules for this 
expansion are as follows, 

• each vertex contributes to a factor ^ and is equipped with an index i taking q values, because of the 
q independent fields in the interaction; 

• each propagator contributes a factor of 1 + v if it connects vertices with the same index and a factor 
of 1 otherwise, 

• sum over all the indices. 

Up to a factor, this is nothing but the partition function Z(f3, J, 7) for the g-state Potts model defined on 
the graph 7. Recall that the latter is a lattice model that is defined on an arbitrary graph 7 by assigning 
spins a v taking q values to the vertices of the diagram. The energy of a configuration of spins a — (a v ) is 
defined as a sum over all the edges of the diagram 

H(a) = -J <k„,<, ( 167 ) 

edges e 

with v and v' the vertices linked by e and J > a constant that favors spin alignment. Its partition function 
is defined as 

Z( / 3,J, 7 ) = ^e-' 3 ^ (168) 

(T 

and it follows from the previous discussion that it is proprtional to the Tutte polynomial 

Z(P, J, 7) = uv v - F Tutto (u + 1, v + 1, 7) (169) 

with v = c~P J — 1. 



4.4 Iterative solution of Polchinski's equation 



Let us now come back to Polchinski's equation (|49| that describes the scale evolution of the effective action 
Sa . Recall that <Sa is obtained from a bare action S\ valid at a very high energy scale by integrating Fourier 
modes between A and A . A large class of propagators that implement this integration and preserve both 
translational and rotational symmetry can be constructed in momentum space as 

p 2 \ ( p 2 \\ S(p+ q) 



A AM ( P .q) =ui Al j-x^;i x (i7o) 
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The cut-off function \ a smooth positive decreasing function with values close to 1 on [0, 1] and close to 
on [1, +oo[. A typical example is x( s ) = e_s that yields a propagator 



A 2 



A AM ( P ,q) = 5(p+q) / i dae~ a P , (171) 

that leads to the convenient expression of Feynman diagrams based on Symanzik's polynomials for 4 theory 
(see [5] and the remark at the end of this section) . Note that we have not included the mass term in the 
propagator to simplify the forthcoming dimensional analysis. Here we treat the mass term as an interaction, 
which does not modify the UV behavior of the theory. The derivative is independent of Ao 

A^^(p,g) = ^x'(^)x<5( P +g). (172) 

and implements the integration over an infinitesimal shell. With such a propagator, Polchinski's equation 



can be cast in the general form (1081 with the non linear operator 



f , , k dA hMl ( S 2 S 5S SS \ 

P{A,S) = / dpdqA A ' ° (p,q) = = = — (173) 

J dk \&4>{p)&m 5ct>(p) 5ct>(q) J 

In high energy physics with c = h = 1, all quantities are measured in terms of masses so that if we 
scale all masses by a factor e s , the action transforms as S — > S' = e sV S with S'[e styD / 2+1 ^ (j)'] = S[4>} and 
4>'(p) = 4>(e s p) in a space-time of dimension D. Then, it is easy to check that 

/3(e s A, e sV S) = e sV p(A, S), (174) 

so that all the results of section |4.1| apply. In particular, a translation invariant term of the form 

J dpi--- dp N S(pi H \-p N ) T(pi, . . . ,p N ) 4>{p\) ■ ■ ■ 4>{pn) (175) 

where V is homogenous of degree n, scales with a factor e sd with 

d= D - N -l\ -n. (176) 

Therefore, we recover the fact that the mass term is always relevant and that for ^> 4 theory in dimension 4 
the coupling constant is marginal and all terms with a high number of fields and/or high order derivative 
interactions are irrelevant. 

Accordingly, Polchinski's equation can be written in integral form, 

,A 

Sa = S Ao + / dA f3(A ,S\>) (177) 

J Ac 

with A = and 

Then, its iterative solution is expanded over ordered Feynman diagrams. The latter are Feynman diagrams 
with a hierarchy of internal lines defined by drawing boxes on the diagram such that the boxes are either 
disjoint or nested and such that a box differ from the next one it is included in by an internal line. Alterna- 
tively, the order on the internal lines can be pictured by a rooted tree t drawn on the set of internal lines of 
the diagram. The contribution A^'^^Sq) to the effective action 5a is computed as follows, 

• associate an intermediate cut-off Aj to any internal line; 

• compute the value of the Feynman diagram using the background field Feynman rules with vertices 
derived from Sq and propagators A Ai ; 
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• integrate the variables Aj over the treeplex T 
Accordingly, the low energy effective action reads 



\ A , defined in section 3.5 



4v°A ( S o) 



(179) 



7ord 



th S 7ord the cardinal of the automorphisms group of the diagram that preserve the ordering. 



Ordered diagrams generate a commutative Hopf algebra Ttp d with coproduct similar to that of Hp given 
in (131 ), except that it has to preserve the box structure. This coproduct is very similar to the rooted tree 



coproduct, as it reproduces on the diagram the admissible cuts of the tree as an extraction of disjoint boxes. 
For example, with A' the non trivial part of A: 




O ® O 



Because the effective action could as well be computed in a single step using the propagator A 



A, An 



E 



E 



^7 7ord S 7ord 



(180) 



there is a Hopf algebra morphism from Tip to 7ip d defined by summing over all the orders on the diagram, 

„ord 



7 



7 



S E/ s 

°7 j i5 'V 

' orders on 7 ' 



(181) 



Moreover, there is also a morphism from the algebra 7i °p A to the algebra of rooted trees which associates n 
ordered diagram with its tree 

(182) 



7 ord -» t 



If we pullback the tree factorial by the composition of these two morphisms we get 

s 7 ^ t\ s; rd ' 

' orders on 7 ' 



(183) 



in accordance with the combinatorial interpretation of the tree factorial in ( 106 1 



From a renormalization group viewpoint, the use of ordered diagrams is necessary since a given diagram 
can be decomposed as different ordered ones that behave slighty differently. For instance, the 2 loop con- 
tribution to effective low energy coupling in </> theory in 4 dimensions depicted in figure 10 can be ordered 
in two kind of ways: Either we consider it a 4 point function inserted in a 4 point function or as a 6 point 
function inserted in a 2 point function. In the first case we get a leading contribution proportional to (log A) 2 
whereas in the second case it is proportional log A. 

To implement dimensional analysis at the level of the Hopf algebras TLf and H™ d we introduce decorated 
diagrams with generalized external structures. Let us decompose the space of action functionals as a direct 
sum £ = £' © £" , with £' containing the relevant and marginal terms, which span a finite dimensional space. 
The decorations are labels i £ I for a basis of eigenstates of the scaling operator V in £' , which we write 
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oc (log A) 2 



oc log A 



Figure 10: Two decompositions of the one loop four point function 



as Oi and an extra index for £ " . For instance, for 4 theory in dimension 4, a basis of £' is given by the 
operators 

(184) 

and we can identify the coupling constant, square mass and wave function parameter as the components 
of the action in this basis. These labels arc generalization of the external structures introduced by Conncs 
and Kreimer in [2], the main difference being that we have to treat all the interactions and not only the 
renormalizable ones. 

A decorated Feynman diagram with a generalized external structure (with an order on the internal lines 
or not) is a couple (7, i) where the vertices of 7 are decorated with indices in /. These diagrams generate a 
commutative Hopf algebra Hj? c with coproduct 

A(7,i) = (7,i)<8l + l<8i(7,i)+ Yl (7i, h) • • • (in, Q ® ( ,~ , , 7 ,~ , , , i) , (185) 



, ~f n disj oint C7 
n decorations £1 



(7i,ii)---(7„,i„)' 



with the vertices of the reduced diagram decorated by the generalized external structures of the subdiagrams. 

Using the background field technique, such a diagram is evaluated with vertices obtained by projecting 
onto the eigenspace of the decoration carried by the vertex for the fluctuating field \ an d the final result 
projected onto the eigenspace associated to i for the background field <fi. Then, all the composition identities 



such as ( 136 1 hold since they simply amount to inserting the decomposition £ = £' £ for each action 
functional. 

The unrenormalized low energy effective action S™ is computed by integrating modes with a propagator 
^4a,a starting with a bare action Sa , so that 

Sr = Va A ' Ao (S Ao ), (186) 

where a is the character that takes the value 1 on all diagrams. The renormalization procedure amounts 
to choose S\ in such a way that at the low energy scale A the relevant and marginal parameters are fixed 
to Sr. If we denote by f3 the character that takes the value 1 on the diagram with relevant and marginal 
external structures, then the bare action is determined by 

Sr^ a ' Ao (S Ao ), (187) 



Then, inverting this equation and substituting in (186 1, we get the renormalized low energy effective theory 
as 

Sr - ^-^(Sr). (188) 

In this discussion, we have used the characters only as a means of selecting some of the diagrams. In fact, the 
evaluation of the diagrams themselves leads to characters thanks to the decorations and external structures. 
This is the point of view adopted by Conncs and Kreimer and is definitely more powerful since it allows 
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Figure 11: A two loop diagram with labelled edges 



to use the algebraic structure based on the Hopf algebra. Unfortunately, in our case it leads to a rather 
complicated formalism when we compose two effective action computations since we have to keep track of 
the location of the insertion of diagrams. 

Finally, let us conclude by a simple example of an application of the previously outlined Lie algebraic 
techniques to the computation of the first Symanzik polynomial. In quantum field theory, a Feynman 
diagram 7 with n edges can be evaluated, in dimension D, as 

d n ry Vy(.<*,p) 

R e"W (189) 
[t/»] 2 

where C/ 7 (a) and V 7 (a,p) are polynomials in the Schwinger parameters a = (ai,...,a„) € K™ and the 
momenta p. Explicitly, the first Symanzik polynomial is expressed as 

c/ 7(«)= E IT a - ( 19 °) 

spanning trees 

where a spanning tree is a tree drawn on 7 that touches all the vertices. 

For example, for the diagram depicted in figure |11| the first Symanzik polynomial is 

ELjA(ar) = 0103 + 0104 + 0203 + 0:20:4 + 0304. (191) 

In the Hopf algebra of Feynman diagrams with labelled edges consider the infinitesimal characters such 
that 

<W~) = 1 and 6 loop (O)=0ti (192) 
and and that vanish otherwise. Then, the first Symanzik polynomial can be written as 

U 7 (a) = e St ™° *e 5l °°p (193) 

Indeed, the first factor extracts all the trees and contracts them. Then, the second factor ensures that the 
reduced diagram has a unique vertex, which means that we extract only spanning trees, and assigns to it 
the product of the a over the non contracted edges. Then, 

U 1 (a) = e 5tr - * e 5l °° p * e~ 5tr00 * e St '^ = e^" 5nla °* * e 8 *™ (194) 

where <5 n i oop is the infinitesimal character that takes the value a i on the one loop diagram with n edges 
and vanishes otherwise. This last relation follows from the Lie algebraic identitity, 

' [^trcc- <5loop] ' ' ' I = <5»iloop (195) 



n iterations 



interpreted in terms of the contraction/deletion relations (4.2). 



3G 



Figure 12: Schwinger-Dyson equation in 3 theory 



Thus, U-y(a) can be evaluated by summing over all contraction schemes of the loops, generated by the 
iteration of the coproduct in the computation of (5 "'°°p. For example, for the diagram depicted in figure 

U ■ (a) — - [(ai + a 2 + a 3 )ct4 + (ai + a 2 + ai)a 3 + (ai + a 2 )(a 3 + a 4 )] (196) 

These terms correspond to the following contraction schemes 

{(123, 4), (124, 3), (34, 12)}. (197) 

In general, such an expression for the first Symanzik polynomial may prove useful in the study of the UV 
divergences, since the latter are governed by the integrals over loops. 



4.5 An algebraic structure related to Schwinger-Dyson equations 

Schwinger-Dyson equations are a set of functional differential equations playing a role analogous to the 
equation of motion in QFT (see [26] for a general discussion). When solved order by order, they allow 
to recover the perturbative expansion in terms of couplings constants. More precisely, they are functional 
differential equations for the derivative of the generating function W[J] of the connected Feynman diagrams. 
W[J] is expressed as a functional of the source J as 

W[J] = h\ogZ[J] (198) 

with Z[J) defined as 

f -SU] + .I -<t> 

Z[J]=M [D0] e—* (199) 





' S ' 











and the normalization constant defined such that Z[0] — 0. Note that we have introduced h so that we can 
measure deviation from the classical theory obtained at h = 0. 

One of the structural properties of the path integral is its invariance under change of variables in the 
space of fields over which we integrate. In the simplest case of an infinitesimal translation <j)(x) — > <p(x) + e(x), 
the invariance of Z implies 

(200) 

Replacing Z [J] by it yields an equation for whose complexity increases with the powers of the 
monomials entering into S. In the simplest case of the 4> 3 interaction, 

S[<j>] = J d D x (-A + m 2 ) + |^ , (201) 

the Schwinger-Dyson equation reads 



d D y J K(x, y)J(y) + +| K(x, y) (J^yj + y K(x, y)- 



(202) 



5J{x) J a \ v ,aj x *' 2 y, »'\5J(y)J 2 v ' y > 6J(y)5J(y) j ' 
with K the kernel of the inverse of — A + m 2 . Let us represent the propagator K by a line, J by a cross and 



W by a blob out of which we draw a line for each differentiation with respect to J. Then, eq. (202 1 is neatly 
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summarized in figure [12] These equations generalize to more complicated interactions: if we consider a 4> n 
interaction we have a term for each partition of n — 1 = n% + • • • + rik corresponding to the product of the nf 1 
functional derivatives of W. Of course, the functional differentials are distributions and the equations may 
become singular when the latter are evaluated at coinciding points, leading again to the UV divergences. In 
this section, we shall ignore this problem, for instance by working with the regularized theory, and focus on 
the formal aspects of the Schwinger-Dyson equations. 



The equation ( 202 1 is a fixed point equation for ip — and may be solved using a geometric series over 



rooted trees. Indeed, rewrite (2021 as 

if = tp + Ttrcc [if] + ^"loop [f] ■ (203) 

and expand tp as a sum over rooted trees with black vertices for Ftrcc and white ones for J-l 00 p- Note that 
Ttrcc is bilinear, so that black vertices have one or two outgoing edges whereas white ones have only one 
because T\ OQ -p is linear. As we compose these operators along the trees, we obtain sums of Feynman diagrams 
with a distinguished external line and terminal lines equipped by crosses. If we iterate Ftree alone we only 
obtain tree-level diagrams while J-\ oov generates loops. 

The correspondence between these Feynman diagrams and rooted trees expansions can be made more 
precise if we introduce a hierarchy on the vertices of the diagram, in complete analogy with the ordering of 
the lines in the iterative solution of Polchinski's equation. The ordering of is encoded in a spanning rooted 
tree drawn on the diagram constructed as follows. The first vertex we meet starting from the distinguished 
external leg is defined to be the highest one and is the root of the tree. As we remove this vertex from the 
diagram, the latter may fall into several connected components. If the diagrams is disconnected after the 
removal of the vertex, we color the root in black, otherwise we leave it blank. Next, in each component, we 
choose as distinguished external line any line that was previously connected to the vertex we removed. Then, 
the lines point towards new vertices that are just below the vertex we removed in the hierarchy. No order is 
assumed between the vertices belonging to different connected components. Then, we iterate the operation, 
considering any component as a new diagram with a distinguished vertex till we exhaust all vertices of the 
diagram. 

It is important to note that the trees we obtain have white vertices with at most one outgoing edge and 
black vertices at most two. Besides, we add extra outgoing edges, called black leaves, to the black vertices so 
that they have exactly two outgoing edges. Then, the tree we obtain are such that below any white vertex, 
there is always less white vertices than black leaves. This property is natural from the Schwinger-Dyson 
point of view: The black leaves represent the sources J and the white vertices the differential operator 
•^loop, which vanishes if there are not enough sources. Conversely, given any tree with this property, we 
can reconstruct all the Feynman diagrams that appear in the corresponding power of the Schwinger-Dyson 
equation by relating each white vertex to a black leaf located below. This tree structure drawn on the 
Feynman diagrams may be used to define another Hopf algebra structure by performing admissible cuts on 
the tree. Of course, the correspondence between Feynman diagrams and trees is not one to one: To a tree 
correspond several diagrams and to a diagram several trees. 

For planar diagrams, the ambiguity in the association of a rooted tree to a Feynman diagram can be 
disposed with by deciding to always move counterclockwise around the diagram as we choose the distinguished 
internal lines, starting from the distinguished external line. This way we associate a unique tree to each 
planar diagram. To reconstruct the Feynman diagram from the tree, we link every white vertex (starting 
with the lower ones) which is not a terminal one to the next black leaf above in the hierarchy of vertices, 
with a link always moving clockwise around the diagram. This way we have constructed a bijection between 
planar </> 3 diagrams with a distinguished external leg and planar trees with univalent (white) and bivalent 
(black) vertices fulfilling the additional condition that the number of white vertices located above a given 
white vertex must be no greater than the number of free edges located below. This is illustrated on a simple 
example on figure |13| Let us note that this bijection is very similar to the one given by P. Di Francesco in 
this lecture [19] . 

The previous bijection is also related to the counting of planar diagrams performed by 't Hooft in |27j . 
Indeed, the sum of those trees with weight 1 is a geometric series and solves the functional fixed point 
equation 

F(x) = x + gF 2 (x) + gh F ^ ~ F(0) . (204) 

Note that here we work with planar trees and count as different two trees that differ by an automorphism. 
This is why we do not divide by the symmetry factor in the geometric series. F is the generating function 
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Figure 13: A planar </> 3 diagram and its reconstruction from a tree 



for trivalcnt planar diagrams with a distinguished external leg 

F(x)= J2 a m g m h n xv, (205) 



m,n,p 



with a mjI i jP the number of such diagrams with m vertices, n loops and p external legs (not including the 
distinguished one). Because of the preivous bijection.The first term in the planar Schwinger-Dyson equation 



(2041 generates the black vertices and the second one the white ones. Besides, the substraction of F(0) 
ensures that there is a black leaf related to the white vertex (otherwise F(x) — F(0) vanishes) and the 
division by x removes this leaf since it is connected to the vertex we have added to the Feynman diagram. 
Because of the previous bijection between trees and diagrams, a m .„ iP also equals the number of connected 
planar (j) 3 Feynman diagrams with a distinguished external leg, m vertices, n loops and p external legs in 
addition to the distinguished one. For example, up to order 3 in g, we find 

F(x) = x + g(x 2 + h) + g 2 (2x 3 + 3hx) + g 3 (5x 4 + lOhx 2 + Ah 2 ) + ... (206) 

If we set h = 0, we obtain the Catalan numbers 1, 2, 5, . . . that count the number of planar binary trees. 



5 Conclusion and outlook 

In this talk, we have presented, in a hopefully pedagogical way, some perturbative algebraic methods common 
to non linear equations and effective action computations. It is conveniently summarized by the following 
table. 



rooted trees 


Feynman diagrams 


non linear analysis 


perturbative path integrals 


fixed point equation 


renormalization group equation 


powers of non linear operators X (x) 


background field technique A 1 (S) 


xi = (id-xy\x) = j: t ^r 1 


= log /[I^e-^^-x+^+x] = E 7 


composition 


successive integrations 



The examples presented here are very simple and can also be treated efficiently using other techniques. 
However, we think that Hopf algebraic techniques may prove to be useful when dealing with more complicated 
situations, in particular in the following subjects. 
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• Multiscale renormalization 



In pcrtubative QFT, we recover a finite theory if we choose a cut-off dependent bare coupling constant 
and expand in powers of the low energy renormalized coupling constant. Nevertheless, the expansion 
is not satisfying since renormalons generally leads to a divergent series. For asymptotically free the- 
ories, the divergence of the perturbative perturbative series can cured by expanding not in a single 
renormalized coupling, but in a scale dependent one The corresponding couplings are obtained 
by performing the renormalization not in a single step but by successive effective action computations. 
Therefore, we expect our techniques to be of interest in the construction of the effective expansion in 
the scale dependent coupling. 

• Graph and matroid polynomials 

Graph and matroid polynomials form a venerable subject where convolution products based on Hopf 
algebras have already proved to useful (see [25]). It may be interesting to further develop the methods 
presented in the derivation the universality of the Tutte polynomial and its relation to the g-state 
Potts. In particular, this may useful to investigate properties of the recently proposed multivariate 
generalizations of the Tutte polynomial [2H], or other graph polynomials like the Martin and flow 
polynomials, which seem to have relations with QFT. 

• Spin foam models of quantum gravity 

Loop quantum gravity (see the book |29] for a general overview) is a tentative quantum theory of 
gravity that provides a construction of the Hilbert space of the quantum gravitational field using 
spin networks. Its path integral counterpart makes use of higher dimensional analogues of Feynman 
diagrams termed spin foams. The latter define a Hopf algebra, as uncovered by Markopoulou [3D], 
which may be useful to address the question of the Wilsonian renormalization of group field theory 
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